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This thesis is a study of the dynamics of a periodically 
supported pipe conveying fluid. The work can be broadly divided into 
three sections ; 

(i) Study of the wave propagation and vibration response of 
a periodically supported pipe ccnveying fluid, 

(ii) Study of the parametric instabilities of such pipes, 

(iii) Investigation of the effects of dynamic absoiber on the 
parametric instabilities of such pipesi 

In the first part of the work, free harmonic wave propagation 
in a periodically supported infinite pipe has been studied, treating 
the pipe as a beanwtype structure. The velocity and the presaire of 
the flowing fluid have been assumed to be constant ihrou^out. We show 



that there erist alternate fre(jien(^ bands of attenuation and propagation 
of free harmonic waves. Due to the presence of the Ooriolis acceleration 
of the fluid, the waves travel in the positive and the negative directions 
with different velocities, numerical results have been computed to show 
the effects of different parameters like the rotational constraint at 
the sujports, the fluid pressure and the fluid velocity, 

natural frequencies of a finite, periodically supported pipe are 
then determined using the ’wave approach'. Due to the difference in 
the velocities of the waves travelling in opposite directions, no 
classical normal modes (in the sense of standing modes) exist. Hie 
effect of introducing identical rotational springs at each support has 
also been studied. The propagation constants have been used* to study 
the divergence instabilities of these pipes. 

The response of a perriodically supported infinite pipe to 
convected harmonic pressure fields has been determined. The resonant- 
type response oocurs at ’coincident frequency' i/diere the convection 
velocity is equal to the phase velocity of the free wave, iHie value 
of the coincident frequency depends on whether the convected pressure 
field travels along the direction of the fluid flow or in the opposite 
direction. The effects of introducing damping in the pipie ^d at 
the supports have also been studied, 

5Ehe second part of the work is corKjerned wi'tti the case 
the velocity of the flowing fluid is having a hsnmonic fluctuation over 
and above the mean constant value. Hiis introduces 'parametric 
instabilities’ in the pipe. Two types of parametric instabilities may 



be distinguidied, the so-called primaiy and seccaadary instabilities, 

ilrst, the method used in previous works for sin^e span pipes have 

been eact ended. Mode shape ajproxLmation is necessary in this method, 

of 

Normally, the approximate modes are taken as those^the corresponding 
beam. Hence, to use this method, mode shapes of a periodically supported 
beam must be known in advance. Moreover, the cciaprtational effort 
increases with increasing number of spans. 

To overc<siie the first pxrobloa, we present a method which does 
not require any mode shape approsamation. This gives a set of coupled 
differential equatiocis, ijy applying boundary conditions, a frequency 
equation is chtained in the form of a determinant. The zeros of this 
determinent gives the boundaries of the unstable regions. The size of 
this determinant, however, increases drastically with the incarease in 
the number of spans of the pipe. 

Another method using the wave approach has also been propwsed 
to detannine the regions of parametric instabilities. Numerical results 
have been presented for the primary and the secondary instability regions 
of a two span pxLpe, An the three methods have been used for the purpose 
of ecmpariscxa and cross-checking. 

The effect of various piarameters like the Coriolis term, the 
flToid velocity and the fluid pressure on the regions of instabilities 
have been sttadied in detail. The effect of damping has also been 


investigated, 



The last jfaase of the work considers the prdblaa of controlling 
the regions of paiametric instabilities by means of dynamic absorbers. 
Eirst^ single span pipes with pinned-pinned, clamped-clamped, and 
clamped-pinned ends have been studied. The effects of various parameters 
such as damping ratio, tuning ratio, and mass' ratio of the absorber have 
been studied numerically, Uo att^pt, however, has been made to optimse 
these jarameters for best performance of the absorbers. Analysis has 
also been presented for the case when an undamped absorber is attached 
at the mid-point of eadi bay of a periodically supported finite pipe* 



CHAPTER 1 


IMEROroCTIOR 

1.1 Inbroduction 

Pipe-lines conveying fltiids at hi^ pressure and velocity are 
encountered in various fields of engineeiing. To name a few, such 
pipe-lines esxist in different process industries, high pressure boilers, 
fuel-lines of air-crafts and missiles, and in nuclear reactors. 

A thorou^ understanding of the dynamics of elastic’ pipes conv^ing 
fluid is recessaiy for the proper design of these pipe-lines. An 
improper design of the pipe-lines may result in leakage and unsatis— 
factoiy performance ■ of the whole system. In severe cases, excessive 
vibration of the pipe-lines may even cause fatigue cracks and eventual 

if 

failure of the pipes. 

In the past two decades, considerable amount of research woik 
has been done in the field of flow induced vibrations. The exact 
analysis of the problem is complicated by the fact that the fluid flow- 
field and the motion of the pipe are affected by each other. In 
theoretical models for pipes with internal flow, fluid flow is usually 
assumed to be undisturbed by the motion of the pipe. Vibrations of 
single span pipes with different boundary conditions have already been 
analysed by many researchers. Some of the woiks analysed the free 
oscillation of a pipe with the fluid flowing at a constant velocity. 
Quite a few studies concentrated on liie parametrfc instabilities of 



a pipe for situations viien the fluid velocity has harmonic fluctuations 
over and above a constant mean value. A historical account of all 
these works can be fcund in reference [ 48] . The effects of bends 
in a pipe have also drawn the attention of some workers [ 24] . 

With slight modifications, the theory developed for the above 
mentioned problons can also be applied to various other problems. 

These include 

(i) the dynamics of a stiff, moving wire or belt [15,52] , 

(ii) the dynamics of amoving chain or string [1,29,58] , and 

(iii) the dynamics of an empty pipe submerged in a fluid flowing 

parallel to the axis of the pipe [ 60 ] , 

With sticng motivation provided by wide applications, mudi 
progress has been made in the study of dynamics of a single span pii>e 
conveying fluid. Very often, long pipe-lines on miltiple supports 
are used to transmit hi^ pressure fluids. However, very little has 
been reported on the vibration problems of sudi multi-supported pipes 
[ 55 ] • The present thesis is an attempt to fill this gap aM deals with 
the dynamics of an elastic pipe supported identically at, regular intervals 
Iterformance of a dyrEuaic absorber in controlling the parametric 
instabilities of such pipes has also been investigated. The objective 
and scope of the present work is outlined in detail in a later section. 



1.2 


Review of Previous Work 



Since the present work deals with the dynamics of pipes conveying 
fluids and analysis of periodic structures, a brief review of the 
available literature in each of these fields is given in the following 
sections, 

1,2,1 lynamics of Hpes Conveying iluids 

Interest in the study of dynamics of elastic pipes conveying 
fluids was activated by Ashley and Haviland [ 2 ] in connection with 
the vibration problems of Trans-Arabian pipe-line. They used an 
approximate power series solution. The equation used to describe 
the motion of the pipe did not adequately account for the governing 
inertia forces, as later pointed out by Reodos’yev [ 15 ] . They 
studied the flow induced damping ina simply supported pipe and found 
it to increase with the flow rate. They further studied the wave - 
type solutions and found that the effect of flow is to expedite the 
waves travelling in the flow direction and to retard the waves 
travelling in the opposite direction. 

Housner [ 23 ] studied the same problem by developing the 
equation of motion using Hamilton's principle. He found that at 
sufficiently higji flow velocities the pipe may buckle, essentially 
111® a column subjected to axial loading, A more general study by 
Uioidson [4o] led to the same conclusions, later on, expeiinaental 
study conducted by Itodds and Runyan [ 14 ] also confirmed the conclusions 
about the stability of the pipe. The experimental results were in 
good agreemert with those obtained by Housner' s analysis [ 23 ] , 



Long [ 28 ] also presented expermeEtal and theoretical studies 
of the transverse vibration of pipe transporting fluid. He used the 
equation of motion as derived by Housner [ 23 ] and a power series 
soluCion was used. His solution was applicable to relatively small flow 
velocities. !I!he pipes with different end conditions were also 
considered. He showed that in contrast to those of simply supported 
pipes, the motion of cantilevered pipes are damped by internal flow 
in a particiilar range of flow velocities. The conclusion of Ashl^ 
and Haviland [ 2 ] that the flow of fluid damps the motion of the 
simply supported pipe was due to improper selection of lesser number of 
terms in the power series solution, 

Hande l ma n 19 ] studied the nature of eigen values from the 
structure of the governing differential equation. No specific solution, 
however, was reported, Movchan [38] used Liapunov’s direct method 
to determine the condition of stability for a simply supported pipe 
conveying fluid, 

Bengamin [3,4] studied the dynamics of articulated pipes 
consisting of rigid tubes connected by flexible goints. He showed 
that a cantilevered ^stem of articulated pipe is subjected to oscillatory 
instabilily. He also found that the dynamical problem is independent 
of the fluid friction. Later Gregory and Eaidoussis [ 17,18] confirmed, 
both theoretically and exi>erimentally, that at sufficiently hi^ flow 
velocities, cantilSvered pipes are subjected to oscillatory instabilities, 

Benjamin [ 3 ] , using his articulated pipe representation, also 
conclnded that bucklic^ instability is possible in the case of a vertical 
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cantilevered system, where gravity is operative. On the other hand, 
Eaidoussis [44] found that vertical continuous flexble pipes never 
experience buckling instability. This contradiction was later 
cleared by Paidoussis and Deksnis [45] . They showed that the 

dynamic behaviour of articixLated and continuously flexible pipes are 
not strictly analogous. The stability of tubular cantilevered pipes 
has also been studied in detail by several workers [ 21,22, 41 ] . 

Maguleswaran and William [39] studied, both theoretically 
and experimentally, the effect of the fluid pressure on the dynamics 
of the pipe. They showed that pipes with both ends supported may 
buckle even at small flow velocities by the action of the fluid 
pressure, Stein and Torbiner [ 57 ] studied the dynamics of an 
infinitely long pipe resting on continuous elastic support. They 
studied the effect of the foundation modulus, fluid velocity, and 
fluid pressvire on various characteristics sudi as the dynamic stability, 
the frequency response, and the wave propagation. 

Srinivasan and Kamath [55] used Eitz method to study the 
vibration of a cantilevered pipe. They calculated the logarithmic 
decrement for the motion of the pipe and also derived the stability 
conditions. The same method has also been used for pipes with other 
conditions [ 56] . 

Thurman and Monte [ 59 ] presented a non-linear analysis for a 
pipe with simply supported ends. They used the perturhation technique 
and found that the importance of the nonlinear term increases with the 



6 


flow velocity. Thus, the applicability of the linear theoiy becomes 
restricted as the flow velocity increases. 

Chen [ 9 ] studied the stability of a pipe with upstream end 
clamped and downstream end constrained by a linear spring, so that 
boiADdaiy conditions are inteimediate between clamped-free and clamped- 
pinned cases. He showed that, in geneiai, both buckling and 
oscillatoiy instabilities are possible depending on the spring constant, 

Ihidoussis and Issid [48 ] derived the most general equation 
for the motion of a pipe. They included the effect of axial contra- 
ction of the pipe. They also showed that with both ends supported, 
the pipe is subjected not only to divergence (buckling) instability 
but may also undergo coupled mode flutter if the flow velocity is 
sufficiently hi^. On the other hand, a cantilevered pipe is subjected 
only to flutter of the single degree of freedom. They also discussed 
the existence of the flutter in conservative system and found it to be 
associated with the gyroscopic (Co2riolis) forces. It was shown that 
flutter in sudi systems can only be of the coupled-mode kind. Some 
design criteria were proposed for maximum flow velocity from the point 
of view of stability. 

In all the above works, the pipe was treated as a beam-type 
structure. Paidoussis and Denise [46,4?] studied the dynamics of 
a thin cylindrical pipe conveying fluid, treating it as a shell-type 
structure. The fluid forces were obtained losing potential flow 
theory. They analysed both cantilevered and clamped-clamped pipes. 
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They fouM that ia addition, to the instabilities in beam modes 
» (corresponding to those fouM previously by treating the pipe as a 
beam-type structure) , instabilities in shell modes are also possible.' 

They also showed that a cl amped-clamped pipe experiences coupled-jnode 
flutter in addition to the divergence instability. The experimeiifcal 
results were in confinnation with the theoiy. Weaver and Unrw [ 6l J 
obtained similar results in case of a pinned-pinned pipe, 

Chen and Roserberg [l2 ] studied the free vibration of cylin- 
drical i^ells conveying fluids. The shell motion was described by 
Hugge's equations and the hydrodynamic forces were derived fi*om the 
linearised potential flow theoiy. In low frequency range, th^ 

developed a modified water-hammer theory and an approximate bouMing 

« ■ 

frequency equation including the effect of fluid flow. 

Liu and Monte [ 27 ] conducted experiments to compare the different 
existing theories. Three velocity regions, viz., low, moderate and 
hi^, were considered separately. Good agreement between theory and 
experiments was obtained in the low velocity region. However, 
discrepancy increased with hi^er fluid velocity. The discrepancies 
were attributed to the fact that the initial configuration of the pipe 
is not strai^t and it also depends on the flow velocity. 

In all the above Judies, the velocity of the fluid flowing 
throu^ the pipe was assumed to be constant. Ghen flO ] studied 
the problem with the velocity of the flowing fluid having a harmonic 
fluctuation over and above the mean constant value. He found that 
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both parametric and combi nati oral resonances are possible. He deter- 
mined the instability regions for a simply supported pipe. 

Ginsberg [ 16 ] also studied the problem of the stability of a 
simply supported pipe with a pulsatile flow. He used Bolotin’s 
method to determine the parametric instabilities. It was found that 
damping, in geraral, has greater effect on the regions of instability 
associated with the hi^er natural frequencies. Similar effects of 
damping was also observed on the hi^er order regions of instability. 

Baidoussis and Issid [ 48 ] determined the regions of parametric 
instability for pipes with various end conditions. They showed that 
cantilevered pipes are also subjected to parametric instabilities, 
provided that the flow velocity is not too small. later Ihidotissis 
and Sundararajan [49 ] determined the regions of conibinational 
resonances for a pipe, ELaut and Huseyin [ 43 ] determined the 
instability regions of a pipe with an axial load. 

1,2.2 Analysis of Periodic Structures 

I^namlcs of periodic structures has been studied by lin [25} 
throu^ calculation of the normal modes and then analysing the forced 
motion in each of the significant modes. Transfer matrix method has 
also been applied for studying the vibration of beams on many supports [26 J, 

Wave motion in periodic systems in the field of crystals, 
electrical circuits, etc., has been studied for nearly 300 years, as 
pointed out by BrillouLn fs] in his classical work. It was orfly 
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recently that wave motion in engineering periodic structures (consisting 
of beams, plates, etc.) has been studied. Heckl [ 20 ] considered a 
system of beams coupled together to form a regular grillage and demons- 
trated that flexural motion in periodic beam type structures can 
propagate in some frequency bands and not in others. 

Mead [30 J studied the free wave propagation in infinite beams 
resting on periodic supports. He also studied their possible intera- 
ction with the acoustic waves, later, he presented a general theory 
for periodic structures with multiple couplings [32] . 

Sen Gupta [5l] developed a graphical method, using wave approach, 
to determine the natxrral frequencies of a periodic becjm. Mead [31 ] 
used the wave approach to determine the response of periodic structures 
wb-ea. subjected to convected harmonic loading. He showed that 
’ coincidence’ occurs at a frequency when the convection velocity of the 
pressure field is equal to the phase velocity of the free waves 
propagating throu^ the structure. Natural frequencies and mode shapes 
of mono-coupled and multi-coupled systems have also been studied in 
detail [ 33,34 ] . 

In all the works mentioned above, periodic element was taken to 
be uniform. ThuSj exact harmonic solution was obtained in closed 
form. But when the periodic element is noa-uniforfe, one cannot find 
the closed form solution. So, attempts have also been made to 
develop approximate methods. Mead and Bijara [36] developed a space- 
harmonic series solution for such problems. Mead and Jfellik [ 35 ] 
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presented an 'assumed mode' method to de-ceimine the response of infinite 
periodic beams va'ch non-uniform elements. Bao and Mallifc [ 5o] extended 
this approximate method to calculate the response of finite periodic 
beams, Finite element techniques have also been used for similar 
problems [ 42 ] . 

Recently, Chen [ 11 ] used the wave approach to determine the 
natural frequencies of a periodic rov/ of circular cylinders surrounded 
by a liquid, 

1,5 Objective and Scope of the Present Work 

It is seen from the available literature, that much progress has 
been made in the study of the dynamics and the stability of sin^e span 
pipes conveying fluids. In the present thesis the dynamics of a 
periodically supported pipe conv^ing fluid has been studied, !Ilie 
woife can be broadly divided into three sections: 

(i) Study of the wave propagation and' vibration response of a 
periodically supported pipe conveying fluid. 

(ii) Study of the parametric instabilities of such pipes, 

(iii) Investigation of the effect of dynamic absorbers onthe 
parametric instabilities of these pipes. 

In the first part of the work, harmonic wave propagation in a 
pexic^ically supported infinite pipe is studied, treating the pipe 
as a beam-type structure. It is assumed that the supports are trans- 
versely rigid and they produce no axial force, Ihe axial contraction 
of the pipe has been assumed to be negligible, The velocity and the 
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pressure of the fluid are taken to be constant throu^out. 

Natural frequencies of a finite, periodically supported, pipe 
are then deteimned using a grajhieal method. In this approada the 
computational effort is independent of the number of spans of the pipe, 
Iho case with identical rotational stiffness at each support has also 
been investigated, lUrtheimore, the propagation constants are used 
to study the diyergence instabilities of such pipes. This is 
followed by the response analysis of a periodically supported infinite 

sub.iected 

pipe^oaconvected harmonic pressure field. The effect of damping in 
the pipe and at the supports has also been investigated. 

The second part of the work considers the velocity of the flowing 
fluid as having harmonic fluctuations superimposed on a constant mean 
value. This introduces ’ parametric instabilities* in the pipe, fhou^ 
the pipe experiences combinational instabilities as well, in the 
present work no attempt has been made to compute these regions of 
instabilities. Two types of parametric instabilities may be distJ-JV 
guished, namely, the primary and the secondary instabilities. Three 
different methods have been proposed to determine the regions of these 
instabilities* 

In the first method, Bolotin’s concept has been used conjun- 
ction with mode shape approximation.- The mode shapes used to 
approxiiaate the displacement of the pipe are taken as those of a 
periodically supported beam. Thus, to use this nethod, the mode 
dbapes of a beam with same number of spans as in the pipe under consd-^ 
deration must be known in advance. 
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In the secocd method, we use Bolotin's concept directly without 
recourse to any mode shape approximation. This gives a set of coupled ‘ 
differential equations. By applying appropriate boundary conditions, 
a frequency equation is obtained in the form of a determmnt.. !I5ie 

size of this determinant^ ! 

however, increases drastically with ah increase in the number of spaiE 
in the pipe, 

Ihe third method uses the wave approach concept developed in the 
first part of the wodc. Instability regions are determined by using 
a graphical procedure similar to the one loentioned earlier. In this 

method, the computational effort is independent of the- mnher of spans 
in the pipe. 

OSie effect of neglecting the mass ratio parameter, g , has been 
studied. We also discuss the effect of damping on the regions of 
instabilities. 

The last phase of the work considers the problem of controlling 
the regions of parametric instabilities by attaching viscously damped 
dynamic absorber to a pipe, lirst, single spjan pipes wdth pinr®d». 
pinned, clamp5ed*-clampied, and clamped-pinned ends have been studied. 

The effects of various parameters such as damping ratio, tuning ratio, 
and mass ratio of the absorber have been studied numerically. Ko 
attempt, however has been made to optimise these parameters for best 
performance of the absorber# Analysis has also been preserved for 
the case when an undamped absorber is attached at the mid-. point 
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of each bay of a periodically supported finite pipe, Por this case 
the instability regions are detemined via the wave approach alone. 

The nomerical results have always been presented in non-dimensional 
forms for various values of the significanb parameters characterizing 
the problems. 



GHAHSR 2 


WAVE PROPAGATIOH AND VIBRATION RESPONSE OE A PERIODIGALIir 
SUPPORTED PIPE COHVEPING EDUID 


2.1 Introduction 

A periodic structure consists of a raimber of identical members 
connected together in identical manners to fotm the whole structure. 
Examples of such structures are, a tall apartment block having 
identical storeys, an aeroplane fuselage structure consisting of 
identical cumred panels reinforced at regular intervals by an orthogonal 
set of identical stiffeners, etc. The identical members forming the 
whole structuire are known as ’ jjeriodic elements'. 

The response of periodic beam-type structu3?es to the convected 
pressure fields has been analysed in the jjast by calculating the 
principal modes of the system, and then studying the forced motion in 
each of the significant modes [ 25 ] . Eor such problems, the 
transfer matrix technique has also been used in which the response is 
obtained directly without recourse to the principal mode calculation f 26 ] , 

When the structure is very long ard contains mary elements, a 
large number of modes must be included in the modal analysis making 
the process extremely unwieldy. This difficulty is avoided, of course, 
in the trarefer matrix approach. However, no clear physical pacture 
emerges out of this approach either. 
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It has heea shown [30] that the steady state vibration of 
periodic structures under the influence of a haimonic source can be 
better understood in terms of a special kind of 'wave motions’. 

These waves are not of simple spatially-sinusoidal form as the supports 
cause reflections and consequent ’near-field' effects. Vibrations of 
such structures to arbitrary excitations can be expressed in terms 
of these wave motions. In particiilar, the response to convected 
harmonic pressirre fields may be readily expressed in simple closed 
form using this wave approach [ 31 ] . 

In this chapter, harmonic wave propagation in a periodically 
supported pipe conveying fluid has been studied, treating the pipe 
as a beam-type structure. The velocity and pressure of the fluid 
have been assumed to be constant througjaout the pipe. When the pipe 
vibrates laterally, the fluid flow field is affected. Conversely, 
the fluid interacts with the wall of the pipe charging its vibration 
characteristics. In the present study, the effect of the fluid 
flow on vibration characteristics of the pipe has been taken into 
account. The effect of vibaration of the pipe on the fluid flow 
characteristics, however, has been assua^d to be negligible, 

Eirst, a brief discussion on the fundamentals of wave approach, 
developed by Mead [30] for a periodically supported beam has been 
presented. Then the receptance method is used to determine the 
propagation constants for the case of a periodically supported pipe 

if : . 

for various values of pressure and velocity of the fluid. Next, the 
propagation constant versus frequency curves are used to determine the 
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natural frequencies of a periodically supported finite pipe. It lias 
also been shown how these- curves can be used for determining the 
divergence instability of such pipe-lines. At the end, the response 
of such pipe-lines to a convected harmonic pressure field has been 
obtained using this wave approach. The effect of damping, present 
in the pipe and in the supports, on the response has also been 
investigated. 

2,2 iUndamentals of Wave Approach 

Let us consider an infinite pipe resting on equi-spaced, identical 

and transversely rigid supports. The pipe is excited at a point by a 

jL03"t / \ 

harmonic force p e (Slgure 2.1), where p and lo are the ampli- 
tude and frequency of the harmonic force respectively, and t denotes 
the time. A form of harmonic wave motion is set up in the whole 
structure, propagating outwards from the point of excitation, IXie to 
the reflections at the supports and the consequent 'near-field' effects, 
this wave motion is not of a simple spatially-sinusoidal form. However, 
on examining the motion at corresponding points in any pair of adjacent 
bays excluding the excited one (for example points H and G- in 
Hgure 2.l), it is found that the phase difference, and the 

amplitude ratio, e , between the motions of such points are same 
for any pair of adjacent bays throu^out the pipe. Thus, in particular, 
represents the phase difference in the motion between two consecutive 
supports, and represents the decay rate over the same distance, 

P is termed as the attenuation constant and p . is known as the phase 
r 1 
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constant. The propagation constant, y , is defined as i yj[_* 

The p 3 ?opagation constant with negative value of y^ corresponds to 
the wave that is travelling from left to ri^t. likewise, propagation 
constant with positive value of y^ corresponds to the wave that is 
travelling from ri^t to left. 

Outside the excited , the wave motion is entirely free and 
is governed only by the free wave equation subject to the boundary 
conditions Imposed by the supports. The propagation constant is 

determined from this wave equation and the appropriate boundary conditions. 
!Eherefore, a propagation constant is a function only of the frequency 
and the pipe/support jfaysical characteristics and is independent of the 
nature of the loading in the loaded bay. Associated with the propa- 
gation constant is a unique mode of transverse vibration of a free bay 
of the pipe. All free bays vibrate in identical modes differing only 
in ihase and amplitude. 

2,2,1 Determination of Propagation Constants 

ELgure 2,1 shows an infinite pipe resting on identical supports 
at a distance S, apart; k^ denotes the rotational stiffness provided 
by each support. Considering any span, say AB in iigure 2,1, one 
can write for free harmonic waves in such periodic structure [ 30 ] 

% = , (2.1a) 

\ , (2,1b) 

■sdiere 6 and M denoteii. the slope and the non-dimensional bending 
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moment, respectively, at sections indicated by the subscripts A 
ani B, Hie non-dime nsi oral bending moment, M , is defined as 
MJJ/EI where M is the bending moment and El is the flexural 
stiffness of the pipe,* Again, 0^ and 9^ can be expressed in 
terms of the moments and throu^ the receptance irelation- 

ships as follows: 

0^ = + 6^^ i 

B BA A BB T 

e = B,, i + 

A AA A AB 3 

#iere is the slope at the section X due to a unit harmonic 

moment at the section Y. Prom equations (2„l) and (2.2), by 
eliminafing ®A» ®B^ \ and Mg , the following equation for the 
propagation constant, p , is obtained: 


(2.2a) 

(2,2b) 




®BB^ " ^BA ” ° 


(2.3) 


Since the receptances B^ are functions of the frequency, 
p is strongly frequency dependent. Obviously, from equation (2,3) j 
for each value of the frequency, two roots for e^ are obtained. 

One of these roots corresponds to the positive going vave and the 
other to the negative going wave. Eeceptances to be used in equation 

(2.3) are derived from__ ^equation of the motion of a pipe conveying 

fluid (see Appendix l). 

If p = p^ + i p^ is a solution of equation (2.3), then 

P = P^ + i t 2Q’'')> n= 0,1,2,... (2.4) 


* Hereafter, all bending mouKcits will refer to the non-dimensional 
values. 
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is also a solution. 02ie pixspagation constants corresponding to 
a = 0 in equation (2.4) are called the primary values of the 
propagation constant. Thus, the wave motion travelling in each 
direction can be thou^t of as wave groups, each having infinite 
components corresponding to integer values of n in equation (2.4). 
Ihe wave components corresponding to n = 0 are called the primry 
waves [ 30 ] . 

As noted in earlier works on beams on periodic supports [30^, 
in the present case also there exist alternate frequency bands of 
attenuation and propagation of free waves. The only difference in 
this case from that of an ordinary beam will be that the values of 
K for the positive and the negative going waves will not be the 
negative of each other. This loss of symmetry is due to the CtoriolLs 
term (due to the fluid flowing in a curved path) in the equation of 
motion of the pipe (see Appendix l). 


Ihen the waves propagate without attenuation, i.e, , with 
= 0, the non-dimensional phase velocities of these waves are 
given by 


and 






i 


for the positive going wave, 
for the negative going wave. 


vhere is the non-dintensional frequency, and and are 

the phase constants for the positive and the negative going waves, 
respectively. 




2.3 natural Prequeneies of a Periodically Supported finite Upe 

Sen Gupta [ 51 ] developed a graphical method based on ^.ve 
approach to determine the natural frequencies of j)eriodically 

supported finite beams. He showed that natural frequencies of an 

# 

H-span beam can be obtained by dividing the range of (in the 
propsagation band) into H equal divisions. Since in the case of 
a pipe, and have different numerical values, a sli^t 

modification of the above mentioned approach is necessary. 

Pigure 2.2 shows an N-sp^an pipe periodically supported on 

transversely rigid supports; the p)eriodic element being a single 

span of the prf.pe with rotational stiffness k at each end (ilgure 

2.3), To determine the natural frequencies, ary dynamic distirrbance 

is considered in terms of two opposite going waves and only propagating 

waves (with p^ = O) need to be considered [ 51 ] . let M_^ and 

M be the moments due to the positive and the negative going waves 

at the extreme left support. If p^ and P^ are phase constants 

for these two waves, then the respective bending moments at the last 
_ iNpt __ iljpT 

support will be e and M e . Since the total moments 

at the end supports must vani^, one obtains 

1 +M =0 (2,5a) 

T* • •• 

iHll iNp. 

M_^e ^ + M_e ^ = 0 ’ (2,5b) 

Suhstitu tion of equation (2.5a) in equation (2,5h) gives 

iNpf iNp“ 
e ^-e ^=0 
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or, 

= 1 2dV1T , 3 = 0,1,2,...,F (2.6) 

Por pidjaaiy values (i.e, , with n = 0 in equation (2.4)) in atigr 
propagation band, the values of 1 between 0 and 2ti. 

He ace admissible values of and to satisfy the boundaiy 

conditions during free vibration, as given by eq'uation (2.6), can be 
TTritten as 

V* = = 237rA, 3 = 0,1,2,... jH (2.7) 

Frequencies within the propagation bands where condition (2.7) 
is fulfilled are the natiiral frequencies of the H-span pipe. 

Thus, once the propagation constants are obtained in the manner 
outlined in section 2.2.1, one can diaw a curve of P* versus 
frequency for each propagatioa band. Ihe value of p* varies 
between 2Tr and 0 in each of these bands. Then, dividing this 
range of p* into H equal divisions, one obtains the natural 
frequencies of an B-span pipe in each propagation band. Thus 'it 
would appear that there are (H+i) natural frequencies in each 
propagation band for an B-span pipe, but this is not so, as explained 
in reference [51 ] , The hipest natural frequency in each propagation 
band, obtained from equation (2,7) happens to be the natural frequemy 
of the individual bays with both ends clamped. Since the ends are 
not clamped in the present case, the hipest natural frequency obtained 
from equation (2.7) cannot be a possible natural frequency. 
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The main advantage of this grajhical approach over the frequency 
determinant method is that the comiutational effort remains independent 
of the number of spans in the pipe. Once the curve of y* versus 
frequency has been obtained, natural frequencies of the pipe with any 
number of spans are determined by dividing the range of p* into 
the same number of equal divisions as the number of spans in the pipe. 

In frequeruy determinant method, however, the computational effort 
increases drastically with the number of spans. This will be apparent 
from Appendix 2 which presents the frequency determinant formulation 
for just two spans. 

It should be noted, however, that at these natural frequermies 
no classical normal modes exist. This is due to the fact that the 
positive and the negative going waves travel at different speeds 
(as implied by different numerical values of and ) and 

hence cannot be supeorfmposed to give rise to standing modes. Thus, 
at these natixral frequencies the modes of the pipe will be complex, 

2,3.1 ITatural Ifequencies of a Unite Pipe on Eqxrl— Spaced 
Supports with Equal Rotational Stiffness 

In section 2.3, jaatural frequencies of a periodically supported 
pipe were calculated using the wave approach. Due to the periodicity 
of the structure, the rotational stiffnesses at the intermediate 
supports were different from those at the ends as indicated in ELgure 
2,2. However, if all the supports are identical, i.e., having equal 
rotational stiffness, minor modification of the wave approach presented 
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in section 2,3 is necessary for determining the natural frequencies. 

Ihis modification is presented in this section. 

Consider an !KUspan pipe resting on equi— spaced and transversely 

rigid supports, each having a rotational stiffness k . If the 

r 

periodic element is taken as a single span of the pipe with rotational 
stiffness k^2 at each end, then extra rotational stiffnesses k^2 
are left out at the extreme ends as shown in ligure 2.4. 

kirst, let us consider ariy one bay of the periodically supported 
pipe neglecting the extra stiffness k _J2. at extreme ends, Hie 
transverse displacemert of this bay, p(^), during free vibration, 
can be written as (see Appendix- 1) 

4 iX C 

yU) = I C e , (2.8) 

np:1 

v4iere the X* s are the roots of the characteristic polynomial given 
by equation (A-1.5), the C^’s are the urknown constants to be 
determined from the boundary conditions, and 5 is the non-dimensional 
co-ordinate along the length of the pipe, 

Por free harmonic waves in the bay under consideration, the 


boundary coioditions are 


o 

II 

o 

(2.9a) 

y(i) =0 , 

(2.9b) 

y*(i’) = y'(o), 

(2.9c) 


where primes denote gifferentation with respect to C , and v 
is the propagation constant. 
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Usiog the above boundary conditions, constants G^, and 
can be determined in terms of C.. Let C =□ G., rt=1,2,3,4, where 

1 7377 

obviously q^= 1. 

Thus, the displacement of ary bay of the pipe, due to the 
positive and the negative going waves can be written separately as 

" + ^ 4 - 

= ^^ I ^ » (2.10a) 

rt=1 

_ 4 _ iA C 

y_(C) = CT I q” e ( 2 . 10 b) 

n=l 

The positive and the negative signs used as subscripts and super- 
scripts refer to the positive and the negative going waves, 
respectively. 

As stated in section 2.3, orrly the propagating waves ( with 

= O) need to be considered for the determination of the nature 

frequencies, let and M be the non-dimensional mcments at 

the extreme left support due to the positive and the negative going 

waves, respectively. The slopes due to these waves at the same 

support are y’^(o) and y*(o), respectively. If and are 

the phase constants for these two waves, then respective bending 

iHpt _ i]jy~ 

moments and slopes at the last support will be e ^ , M__ e ^ 

ilTP^ _ iHp" • 

aM (o) e , - y’ (o) e * Hie total bending moments at 

the end supports will no lorrger be zero, as given by equations (2.5). 

Due to the presence of extra stiffness k^2 at these ends, the 

modified boundary conditions are 
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ty'(o) +y»(o)} , 


__ ilTy^ iffy" < iHyf i]jy7 

M^e +M_e ^=-|{y^(0) e +?(o) e ^} 


(2.11a) 


(2.11b) 


vHaeve is the noii-dimensional rotational stiffness, 

fhe bending moments and M can be written as 


= - y+ (D) + ~ rjio) , 


(2.12a) 


M_=-y"( 0 ) + y^( 0 ) 


(2.12b) 


Substituting equations ( 2 , 10 ) and (2.12) into equations (2, II), 


one gets 


i, ’n i * <=■ ^ (’'n + i “ ° (2- '5^5 


- /.2 


+ 4 iWy" 4 - 

3 . e Yq ^^+Ce ^yq”x ^=0 

1 J. . 'in n "^1 I % A u 

n?:l n=l 


( 2 . 13 b) 


Por non-tri-vial solution of C and C~ 

1 1 

I 4 (4 ^ "r’' J I 4 (’•n * ^ V’-n) 

n=1 n=1 


^®4 i + 2 

* I %>-a 

n =1 


iuy. 4 „ 

e ^ Z q" x 2 

ii “ “ 


a a 

e ^ [ I q* ] [ T q- (X^ + i K X ) ] 

'■ ^ ^ ^n n n ^n n r n-^ 

n=1 n=l 


4 . . 

1 _ r + /, 2 


— e rlt (^ +i*‘‘^)iry q i ] 

L i /^n ^ n r n'^ ^ ^n n 

n?=l n ?=1 . 


(2.14) 
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The left hand side of equation (2. 14) is strongly frequency 

deperdent and in general, it will be complex. Erequencies in each 

propagation band, where equation (2. 14) is satisfied, are the natural 

frequencies of the pipe. For this case, the grajfaical approach 

presented in section 2.3 can no longer be used. However, once the 

propagation constants (and corresponding ^ *s and q ’s) are 

n n 

obtained, natural frequencies for ar^r number of spans in the pipe can 
be determined using equation (2,14). It is interesting to note 
that in this case also the computational effort is independent of 
the number of spans in the pipe. 

Ifetural frequencies of a periodically supported finite pipe 
on rigid supports without rotational stiffness can be determined 
considering = 0 in either of the methods presented in sections 
2.3 or 2.3.1. 

2.3.2 Buckling ( Diverge nc e) Inst ability 


A careiUl examination of equation (A-1.2) reveals that the 

2 

coefficient of the term is equivalent to an axial load on the 


pipe. Hence it is apparent that an increase in the velocity or 
pressure of the fliiid beyond certain critical values will cause static 
buckling of the pipe [ 48 ] . This static buckling is also termed 
as "Divergence Instability", Bolotin [ 6 ] has shown that static 
buckling can be considered as a special ease of dy nami c iretability. 
Buckling is implied when at certain critical values of the parameters, 
one of the mtural frequencies of the ^stem is reduced to zero. 
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Hatiiral frequencies of a peiiodically supported pipe can be 
obtained following the method developed in section 2,3, If, for 
certain combinations of parameters any one of the natural frequencies 
is zero, it will imply the onset of static buckling of the pipe-line 
in the corresponding mode. Hence, combinations of fluid velocity, 
fluid pressure, and rotational stiffness, which make the pipe to 
buckle in a particular mode can easily be determined. 


2.4 Response of an Infinite Periodically Supported Pipe to a 
Convected Harmonic Pressure ELeld 


Mead [31 j studied the response of periodically supported beams 
to a convected harmonic pressure field using the wave approach. He 
showed that coincidence occurs at a frequency when the velocity of 
the free harmonic waves travelling in the beam is equal to the velocity 
of the convected loading. In the present section, the response of a 
periodically supported infinite pipe to a convected harmonic pressure 
field has been studied. 

Consider an infinite periodically supported pipe resting on 
transversely rigid supports and excited by a convected harmonic 
pressure field which exerts a force 


f A.\ - kx) 

pjx,t) = P„ e ^ 


(2.15) 


per unit length of the pipe, !Hae convection velocity, G^, is 
/k and the phase difference between loading at points a distance 
^ apart is e = -k& , 


The equation of motion for transverse 
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vibration of the pipe, between argr two consecutive supports, is (see 
Appendix l) 


El -2^ + (b^ T^+ AJ H + 2 -¥■ 


1 ■ v-f ■ ■ j-j ~p^ 2 ■ ‘ “f ' 3x8t ' '“f “P' ,.2 

OX ^ 3 x at 


= TD e' 


f ‘ 9x3t 
i(ojt - kx) 


( 2 . 16 ) 


On substituting y - Hy , equation (2.16) can be pit in 


the following non-dimensional foim 


+ (u^+y) + i2 6 uQ^ - ^^y = p 

m dr ^ 

p 

v^iere the non-dimensional loading parameter, Pq= — 


(2.17) 


e can be written as e = — fi/ C , where C is the non-dimensional 
convection velocity of the loading^ defined as 

C = ( 0^ (2.18) 

Okie solution of equation (2.1?) is 


4 ix 

y(c) = I C^e + 


- ieC 


J3?=1 


4/2x2 2 

e — (u +y) £ — 2^viQe - Q 


#iere the X Vs are the roots of the polynomial 
n 


(2.19) 


X^ - (uVy) X^ - 2 0unX - = 0 (2.20) 

Bie C^’s in equation (2.19) are to be detemined from the 

boundaiy conditions as discussed below. 

Consider any two adjacent elements of the pipe as ^own in 
Hgure 2.5. ^ virtue of the forced wave property, the slopes and 
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momenfcs at — end r of the element II and at end r— 1 of the 

element I are related as [ 31 ] 






(2.21a) 


M 


r 



(2. 21b) 


■wtoere primes denote differentiation with respect to 5 . Now, the 
displacements at the supports are zero, i.e. , 


= 0 J 

^r-r ° 

lihe bending moments and ^ of equation (2.21b) can be 

written as 


K 


Mp = - y'’(i) - -f 

?(l) , 

(2.23 a) 

f 

II 

1 

+ 

y’ ( 0 ) 

(2.23b) 


Use of equations (2.19), (2.2l), (2.22) and (2.23) yields the 
following system of equations for the coefficients C^'s used in 
equation (2.19). 


I °n ~ “ 4 / 2 s 2 ^ 2 » 

n=1 e -(u +y) e - 2gu fie 


4 iX 

I c e ^ 

li=1 “ 


- ^xe 


4 / 2 V 2 2 » 

e -(u +y)£ - 2guQe 


4 iX 4 

I 0 X e “ = y 0 X 

an a n 

0=1 0=1 


(2.24a). 


(2.24b) 


> 


(2.24c) 
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4 

I 

a=1 


C (x^ - 
n n 


2 


iXj3^) 




ie ^ 

= I 

Et=1 


c (x" 

n n 


K 

+ ~iX ) 
2 n.' 


i K e p 

+ "4 2" ' ^ ~2 ° 2 ^ (2.24d) 

E -(u +Y)e -20u£2e-n 

Oace the have been obtained from equations (2.24a)- 

(2.24d), the response of the pipe to the convected harmonic pressure 
field can be determined from equation (2.19). Hysteretic type 
damping in the pipe and at the supports can easily be incorporated 
by replacing EL inequation (2.1 6) by Sl(i+ ip^) and in 

equation (2.23) by *^^.(1+ where and are the loss 

factors of the pipe and the supports, respectively. It should be 
noted that the introduction of damping will change equations (2.17)} 
(2,19)} (2.20) and (2,24) accordingly. 

2.5 Eesults and Discussions 

2.5.1 Computations iterfomed 


The following computations have been performed using the analyses 
developed in previous sections. All the results a3?e expressed in 
terms of non-dimensional quantities. 

(a) Curves of propagation constants, y (=yp'‘" 'versus frequency, 

n, have been obtained for an uoiamped pipe with different conijinations 
of fluid velocity, u, fluid pressure, y , rotational stiffness, ar^ 

mass ratio parameter, g . Equation (2.3) was used to determine the 
propagation constants. Basse velocities of the prLmaoy waves have 
also been computed. 
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(b) Mimerical values of the first two natural frequencies of a two 

span, pipe, with ~ ^ f have been determined using the method 

discussed in section 2,3. As a check, same values were computed 
using the frequency determinant given in Appendix 2. 

The first two natural frequencies of a two span pipe of the type 
discussed in section 2,3.1 have been calculated using equation (2.14). 
Valuations of these two frequencies with the value of have also 
been shown. 

Various combinations of the first critical values of the fluid 
velocity and the fluid pressure to cause buckling of the pipe have 
been piresenbed. 

(c) The response of an infinite pipe, to convected harmonic pressure 

field with different convection velocities, has been obtained using 
equation (2,19). Bae response is expressed in terms of the non- 
dimensional curvature p (= ) at the mid-point of a bay. The 

effect of variation in the damping of the pipe, and at the supports, 

n , on the level of the response has also been discussed. dSie non- 
dimensional loading parameter, p^, has been taken as unity in all the cases, 

2.5.2 Propagation Constants 

Using equation (2.5), the propagation constant, y , has been 
calculated for different values of the frequency, n . Sigure 2.5 
i^ows the variation of y with Q for u=2., Y =2, 0=0.5 and k^=0, 

When y =0, free waves can travel throu^ the pipe without attenuation 

r ^ • 

and the ranges of over whidr iij, = 0, are called the ’’pEropagation 

It is seen from Tigure 2.5 that there are alternate attenuation 


bands". 





a)3d propagation bands. It is found that the lower bounding frequency 
of the nth propagation bard is identical with the nth natural frequency 
of the single bay periodic element shown in iigure 2.5« Hie upper 
bounding frequency of the nth propagation band is identical to the nth 
natural frequency of the same periodic element v^hen its ends are 
clamped. It is clear from iigure 2.5 that in general, the values of 
and are different, Ihis is due to the Coriolis acceleration 

of the fluid while flowing in a curved path. As the Coriolis term 
is proportional to V (and not to 7 like the centrifugal force term), 
the system loses symmetry along x-direction, Herxje, waves travelling 
in the direction of the fluid (called positive going waves) have 
different velocities then those travelling in the opposite direction, 

The values of the jhase constant, p^, vary in the attenuation bands. 
However, the variation is found to be decreasing with successive 
attenuation zones. 

To study the effect of parameters like k^,Y, u, and 3 on 
the propagation constant, calciXlations were performed for different 
combinations of these parameters. These results, shown in Eigures 
2. 5-2, 8, can be summarised as follows: 

(a) Eigures 2.5 and 2.6 show the propagation constant for different 

values of < , the other parameters being same in both the oases, 
r 

It is seen ft am these grajhs that the effect of increasing the rotational 
stiffness, k , is to shift the start of the propagation bands to hi^er 
frequencies. However, the ends of the propagation bands ramain 
unaffected. As noted earlier, the lower bounding frequencies of each 
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propagation band correspond to the natural frequencies of the single 

bay element of iigure 2.3» Ihat is vdiy, these shift towards higier 

values with increasing value of However, the upper bounding 

frequencies of each propagation band coirespond to the natural 

frequencies of the single bay element of figure 2.3 'vdien its ends are 

clamped and as suda remain inaependent of the value of ic . 

r 

(b) ligures 2;6 and 2.7 show the effect of increasing the fltiid 
pressure, Y , with the other parameters remaining the same. Hie 
curves show that the effect of increasing the fluid pressure, y , is 
to shift both ends of the propagation bands to lower frequencies, 

(c) iigures 2.6 and 2,8 show the effect of increasing the fluid 

velocity, u. The increase in the fluid velocity, like that in the 
fluid pressure, shifts both ends of the propagation bands to lower 
frequencies. However, the fluid velocity has more pronounced effect 
than the fluid pressure. It can be seen from Iigures 2.6 and 2.8 
that increase in the value of u from 2 to 4 shifts the end of the 
first propagation band from 20.6 to 0 = 16.4 approximately. 

On the other hand, for the same increase in the fluid pressure, this 
drift is from 0= 20,6 to = 19,8 only, as seen from iigures 2.6 
and 2,7. 

(d) niie curve of v versus changes very insignificantly with 

changes in 6 in its physical range, i.e., 0 1. S 1.1. Keeping in 

mind that the nonr-dimensional frequency, Q , includes m^, the above 
statement should not, however, be confus,ed with the effect of the fluid 
mass, on the value of the propagation constant. 
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ELgurG 2* 9 shows th© jhasG velocities of the ptimaiy waves 
propagating in the positive and the negative directi ons» 0?o save 

space, the negative of the phase velocity of the negative going wave 
has been plotted* It is seen from the g 3 ?aph that for prdmaiy waves, 
the velocity of the positive goirg wave goes to infinity at a 
frequeccy when = O. Beyond this frequency, the sign of yf 
changes, so there is no more positive going primaiy wave. Both the 
pa?imaiy waves travel in the negative direction and their velocities 
finally attain the. same finite value as indicated in figure 2.9. 

2.5.3 ISIatural Frequencies 

As discussed in section 2.3, to determine the natural frequencies 
of an It-span pipe, the grafh of y* (= |yT - yTj ) versus ft is 
plotted in Figure 2.10, A frequency range covering only the first 
propagation band has been considered. Values of yT and y^ 
have been taken ffom Figure 2,5. Sbr an M-span pipe this range of y* 
is divided into if equal intervals. The first N points of this 
division correspond to the first H natural frequencies of the pipe, 
iirst two natural frequencies of a two span pipe with simply supported 
ends ara indicated in Figure 2,10. These values come out as 6*11 

arri ft *s 12,82, and are in very good agreement with those obtained 
from the frequenegr determinarit, discussed in Appendix 2. Ehe sii^t 
error can be attributed to the difficiilty of reading the values very 
accurately from the graph. It has also been Checked that 
first natural frequency of the sin^e b^ vhen one of the ends is pinned 
and the other clamped. As stated earlier, ft^^ is the first natural 
frequency of a single bay with both ends simply supported. 







To study the effect of rotational stiffness, on natuial 
frequencies of a pipe resting on identical supports at regular 
intervals, first two natural frequencies of a two span pipe have heen 
determined using equation (2.14). Tigure 2.11 shows the effect of 
on the first two natural frequencies of the pipe. As expected, 
it is seen that increase in increases both the natural frequencies. 
However, the difference in their values reduces with increasing value 
of la the limiting case when is infinite, both 

frequencies converge to the first natural frequency of a clamped— 
clamped pipe implying independent vibration of the two bays. 

As discussed earlier in section 2,3.2, if any one of the natural 

frequencies of the system is reduced to zero, it implies the otBet 

of the static buckling in the corresponding mode. Moreover, the 

first natural frequency of a periodically supported pxipe with any 

number of spans, is given by the start of the first propagation band* 

The combinations of the fliiid pressure, Y , and the fluid velocity, 

u, for a given value of the rotational stiffness, are determined 

which make the first propa^tion band to start at zero frequency. 

In other words, the miaiirum values of Y and u necessary for the 

onset of static buckling of such pipie-lines are obtained, !iaiese 

results are shown in Etgure 2.12. It is found that the equation of 

2 

the curves shown in iiguie 2,12 is given by u +Y= corstant, where 

the value of the constant depends upon that of ¥alues of this 

2 2 

constant for k =0 and » are ir and 4v , respectively, 
r 

The analogy between the buckling of a pipe and -Kiat of an ordinary beam 
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can be sou^t by taking axial load on the beam equal to u^+ Y; and 
= 0 corresponds to the pinned-plnned beam and = » to the 
clamped-clamped beam. 

2.5.4 Eesponse to a Gonvected Harmonic Pressure Eield 

Using the analysis presented in section 2.4, the response of a 
peiiodically supported infinite pipe to a convected harmonic pressure 
field is calculated. The response has been expressed in terms of 
the non-dimensional currattire amplitude, p , at the centre of any bay. 

iigure 2,13 shows the response of the pipe for different velocities 
of convection, C, of the pressure field. A frequency range covering 
ohly the first propagation band of the corresponding \indamped pipe 
has been considered. Prom Pi&ire 2,13, it is seen that the response 
attains maxima at certain frequencies. Inspection of equation (2,19) 

I 

shows that this will occur if 

(a) ^ assume very large values. This happens ^iien the phase 
difference, e , between pressures at distance z apart is equal to the 
phase constant, for free wave motion (in absence of damping). 

In other words, resonant -type response occurs when the convection 
velocity of the pressure field equals the phase velocity of one of 
the free wave components. This phenomenon is termed as ’coincidence’, 

(b) The denominator of the second term of the rl^t hand side of 
equation (2, 19) has absolute minimum value ^ich will be zero if 

there is no damping in the pipe. The frequency at vdaich this occurs 

-2 

is the so-called resonant frequency given by n = C . 


Depending 



50 


upon 1:110 value of Cy coincidence may occur af a frequency lower fiian 
the resonant frequency. 

As observed in the case of a beam, Slgure 2.13 also shows that 
an increase in the convection velocity shifts the coincidence frequency 
to hi^er values [ 31 ] . 

It should be noted that for the pipe, the free wave phase velocities 
are different for the jxjsitive and the negative going waves. Hence the 
coincidence frequency depends on whether the convected pressure field 
travels along the direction of fluid flow -(as as sume d ia section 2.4) 
or in the opposite direction. 

ilgure 2,14 ^ows the effect of damping in the pipe (n^) on the 
level of the response. It is seen that the damping suppresses the 
response and peaks are reduced considerably. lampiEg in the pipe is 
effective in reducing both ’coincidence’ and ’resonance’ peaks. 

Ilgure 2. 15 shows the effect of damping in the support (’1^) on the 
response of the pipe. It is seen that only the ’ coincideiKje’ peaks 
are effectively reduced by the damping present in the supports. 

However, the resonant peaks are not affected mudi by . 





CHAITEE 3 


PAMIilEffEIC IHSPABILITIES OE A PERIODiGALIff-SUPiOHrED PIPE 

COFyEyilTCJ ELUID 

3.1 lat reduction 

It is well known that if the velocity of the fluid flowing throii^ 
a pipe is not constant and is having harmonic fluctuations over and above 
a constant mean value, then the pipe experiences parametric instabilities. 
This phenomenon is similar to that of a beam when subjected to a periodic 
axial force [7]. Ihrametric instabilities differ from the more 
familiar instabilities , such as elasto-static buckling and forced vibration 
response, in the following aspects : 

(i) It occurs over a range of parameter space and not at discrete 
points. 

(ii) It can occur in a direction noimal to the excitation. 

(iii) It may occur at frequencies other than the natttral frequencies of 
the system. 

Two types of parametric instabilities may be distinguished, namely 

the primary and the secondary instabilities. For a system with natural 

frequencies given by (n = 1,2,3,.., )f prinaiy instabilities occur 

2 

at fi s= ~ n with K - 1,3,5,..., as the excitation parameter is reduced 

K a 

to zero, where n is the frequency of the periodically varying parameter, 

2 

likewise, secondary instabilities occur at n = ^ with K = 2,4,6,... 
as the excitation parameter is reduced to zero. Instabilities corresponding 



to K — 1 ard 2 are known as the principal primary and the principal 
secondary inscability, respectively. Bolotin [7 ] has shown that of 
all the instabilities, the principal primary instability is the most 
critical and has, therefore, greatest importance. Ihis region of 
instability is also known as the principal region of dynamic instability. 

The instability regions of a single span pipe with different 
boundary conditions have so far been st idled in detail [ 10,16,48] . 

In all these works, comparison functions'^ [j?] were used to represent 
the motion of the pipe. Thereafter, Bolotin's method [ 7 ] was used 
for the determination of the instability regions. In this chapter, we 
discuss the parametric instabilities of a periodically supported pipe. 

The supports are taken as transversely rigid and without ary srotational 
stiffness, Bor the determination of the instability regions, following 
three methods have been attempted : 

(i) The method used in earlier works (mentioned in the above paragraiph) 
has been extended for the present piroblem, 

(ii) Bolotin* s concept is used directly without taking recourse to ary- 
mode shape approximation. In this method, the problem reduces to 
finding the zeros of a determinant. The size of this determinant, 
however, increases drastically with the number of spans in the pipe. 

(iii) A method has been developed using the concept of wave approach 
presented in Chapter 2. In this method, the amount of computation 
is independent of the number of spans in the pipe. 


t sa-cisfying all the bomdaiy conditions (geometric and natural). 
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numerical resitLis have been computed for the piinHiy and the 
secondary instabilities associated with the first two modes of a two 
span pipe, la all cases, only the principal regions have been determined. 
The results obtained by all the three methods are found to be in excellent 
agreement with one another. The effects of various terms, like icass ratio 
parameter, 6 » mean fluid velocity and fluid pressure, on the regioce of 
instabilities have been studied. The effect of damping in the pipe has 
also been investigated. All the instability regions are ^ own as plots 
of the nonr-dimensional frequency, fi , versus the excitation parameter. 

The excitation parameter has been taken as the ratio of the amplitude of 
the harmonically varying component to the mean value of the fluid velocity. 


3,2 Eqmtion ibr Transverse Motion of the Eipe 

Consider an liUspan pipe idenbically supported on transversely rigid 
supports at distance apart. The rotational stiffness at each support 
has been assumed to be zero. The equation for transverse motion of 
the pipe in ary span is (see Appendix l) 


El + (m^ + PfAp) ^ + 2mf ^ + “f H ^ V ° 

3x 


>5 


2 

J 

(3,1) 


Equation (3.l) can be pub in the following non-dimensional form ; 


3 y. 


2 - 

3 y. 


2 - 

3 y 


sy. 


2 - 

3 y. 


1 ? 3 yi ^ ^1 au '’•'1 ° ''I 

i + (u^ + y) — ^ + 2 eu + e 17 ~ 


3 t 35 


3t2 


(3.2) 


ar 

where the non-dimensional displacement of the pipe, y ^ , the non- 

t 
£ 


dimensional time, x = ( ~-) \ , and other symbols have already 
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been defined in Chapfter 2. 

Let the aoii-diffl.ensioaal velocity of the fliiid, u, be of the form 

u = Uq (l + 6 cos fh-) , (3,3) 


where is the constant mean velocity of the fluid, 

<S is the excitation parameter, 
and n is the non-dimensional frequency. 

Substituting equation (3,3) into equation (3,2), one gets 


4- 2- 

3 y. 2 2 3 

— ^ ~~o + 6 cos fix) 


3C 


95 

37 


fii 

3C3t 


- 6 n. fi6 sin fhr 


1 


95 


2- 

3 y-i 
2 

3t 


= 0 


(3.4) 


Since equation (3.4) contains periodic coefficients, it may have 
unbounded solutions even for small value of the excitation parameter, 

5 [7 ] . The unbounded solutions will occur over a range of frequencies. 
These regions of instabilities are of interest to us and can be 
determined as discussed in the following sections. 


3,, 3 DetermLoation of Regions of ifeiametric Instabilities 


Bolotin [ 7 ] has shown that the regions of unbounded solutions, 
of equations of type (3.4), are separated from the regions of boiaaded 
solutions by the periodic solutions with periods T and 2T, where T = 2-n/^i 
(see Aj^jendix 3), In other words, two solutions of identical periods 
bound the 3 ?egion of instability, and two solutions of different periods 
boutd the region of stability. The regions enclosed by solutions having 



period 2T, correspond to ' pxlmry instabilities', vdiereas, 'secondary 
instability’ regions are enclosed by solutions having period T. 

Hires different methods to determine the solutions of equation (3.4) 
with peinods 21 and i.e,, the regions of primaiy and secondary 
instabilities, are discussed below. 


3.4 I'irst Method 


In previous works, Bolotin’s method was used in conjunction with 
mode shape appro 3 d.mation for determining the regions of parametric 
instability of a sia^e span pipe with different boundary conditions 
[ 10,16,48] . In this section, the same method has been used to compute 
the instability regions of the pipe considered in section 3.2. Hie 
transverse displacement of the pipe is approximated by the mode i^apes 
of a periodically supported beam. Hence, to use this method mode shapes 
of an H-span beam must be known in advance. 

Let the displacement of the pipe, y^ in the transverse direction 
be approximated as 


= I Qt^CT), (3.5) 

r 

q^’s are the generalised co-ordinates (functions of tir^), acd 
ar6 the couparison functions satisfying the bouniaiy conditioi^^ 


Foi* the in quest ion^ ^ 's can be taken as the nomB-lised modes 

4 

of an HUspan beam# H^nce, ^ satisfy the equation — — “ 
Substitirfcing equation (3«5) into equation (3«4)> one gets 
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I [ ^ u^(i + 5 cos J2 t) <|)^ + {<{)^'' + (u^ (l + 6 cos Qt)^ + y) 


- 3 Uq fi(S sin CJx $p ]= 0, 


(3.6) 


where primes and dots denote differentiation with respect to C and t , 
respectively. 

Multiplying equation (3.6) by vf and integrating over the nocK 

JM S 

dimensional length of the pipe, JST, one gets 

Q + 2$ Uq(1 + 6 gos nx) [P] Q + { [S ] + (u^ (l + fi cos Or) ^ + y) [R] 


- 3 [^1 ^ Q = 0, 


(3.7) 


where Q is a vector J q 2 i j and the elements of the matrices 


P, R, and S are given by 


1 

^mn R / ‘*’n ’ 

0 


r = J- / A A n gr 

^mn R ' ^m ^n ^ ’ 


xt R 




fi (by im>-ing use of orthogonality of the modes) , 
''n mn ^ *' 


where S = 1 if m = n 
mn 


= 0 if m n 


Thus, Tibe matrices, P, E, and S can be deteimined if the mode shapes 


are known. 
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As stated earlier, the regiocB of primary instability are separated 
by solutions havirg period 21. itor these solutions, vector Q can be 
written, as 

Q = I (-2 +bg. cos , (3.8) 

K=1,3,5^--- 

Tfltiere a^. and b^. are unknown vectors. • 

Substituting equation (3.8) into equation (3.7), one obtains 

u^ ^2 

I f V ^""o “V"^ tR] ) • 

K=1,3,5^--- 

sin (-I KQx) 

2^2 

+ {-(p^ + 2e Uq [P] (P + ([S] + (u^+Y + ) [l(l )b^}cos(^ Kftr) 

+ {-61^0 <5 [^] (i^)5K + ^ “ 2 

+ {6 u^ 5 [P](~) fK ^ 2 ^ 

+ { - eu^ 6 [P] (~) bj^ + u^ 5 [ e] a^+ I e u^ 5n[p] b^} sin Or) 

+ {6Uq 5 [P] (~) ^+%6 [e] bj. - ^ e u^SQ [I] ^ > cos Jit) 

u2 ^2 

+ — -- |— [e] {hg sin (^^ Qx) + bgcos (^|^ Ox) + ^sin (^~ fix) 



COS 
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where 




D 

11 

12; 

t 

13 


D ; 

D 

21 

22 1 

23 


^2 

% 

D, . 

D. 

D 

41 

42 

43 

• 

• 

• 

• 

• 

• 

• 

• 

« 


®14 

' a. : 

1 -1 1 

^24 — 

!b, I 

! ■ i 

^4 •-•** i 

i ! 

{ { ? 

”44 

lb ^ 

; ^ i 

• 

• 

!* I 

1 ■ 

♦ 

i • : 

i • '* 


= 0 , 


I >11 =-r[i]+[s]-(uJ+r*^-u 2 s) [H], 

d .9 = "6 u n [p] , 

2-2 

“l 3 = ’ 

= - 6 6 n [ ? ] , 

Dj, = 6 !! i I ] , ^ ^ 

D22="^[i] + [s] + (u^-t-Y+ °- - + Ugg) Ce 3 , 

®23 5 « [^1 » 

2 2 

“24“ ("^+-^ )C=3 » 

- 


^2 " ~e\«n r^p] * 

22 

1^3 .- 2 £[i] . [s].(u= [b], 

%4 =- 3 Su^fl [P] . 

^41= 8u^ «a [p] , 


(3.10) 
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V = tB], 

®43 " 3Bu^n [P], 

oq2 2 

®44 [I] +[S] + +«+-V) [H], 

and [ I ] is the identity matrix. 

Preqiiencies bounding the instability regions are obtained by 
setting the determinant of the coefficient matrix, [ P] , of equation 
(5«10) equal to zero. Of course, the de'cerminant is of infinite 
order, but it belongs to the class of normal determinants and is 
absolutely convergent [ 7 ] . Hence, the approximate boundaries of 
the instability regions may be obtained by truncating the series in 
equation (3.8) at K = 1. Oliat is, by equating to zero the determinant 
shown circumscribed with da^ed lines in eqi^tion (3.IO), !Ehis is 
known as K = 1 approximation, \shicii cecessarily yields only principal 
region of instability, A better approximation and as well as hi^er 
order instability regions would be obtained if the series in equation 
(3.8) is truncated at hi^er values of E, Of course, the series in 
equation (5.5) must be truncated at adequately hi^ value of r, which 
defines the order of P . Bolotin [7] has shown that in the ease of 

ii 

a beam subjected to a periodic axial force, for values of 6 even upto 
0.3, the regions of prindial primary instability obtained by t = 1 
approximation, are within one percent of those obtained by K = 5 


approximation. 
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low, the regions of secondaiy instability, ishidi are separated 
by solutions of ^riod T, may be obtained by expressing 

S= I + b cos , (3.11) 

K=0,2,4_,'”- 

which substituted into equation (3. 7), yields once again' equation 

(3*9), but with sunmation now being over K = 0,2,4,.... Equatii^ 

coefficients of similar terms as before, one is led to a matrix 

b^^ 

equation* equivalent to eqxjation (3.10), but with a vector i * aM 

i 2 •, 

i^l 

0 

'■ * 

thence to a vanishing determinant which yields the boundaries of the 
secondary instability. Olie dexernananb obtained for K = 2 approxi- 
mation is given by 

®13| 

D23 , (3.12) 

%3 

where ^11 “ [ ^ 1 7+ *“2 )[®1 » 

©12 = ^ e % » 

^15 =u^«[H] , 

^2^2 

D22=[S]‘^ (Uq+y + -\-) [h] > 

D25 " 

3^1= 2 ul 5 [R], 

D 2 = 2auQn [p] i 

1^3 = [s]+ (uj [B] 
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3.5 Secoad Method 

la this section, Bolotin' s concept is directly used without 
recourse to ai^y mode shape approximation. To detettnine the regions 
of primary instability, vihich are bourded by periodic soltcdions with 
period 2T, the displacement, y^, can be expressed as 

= I i^(5)sin(^T) +Yj^(C)cos(^fiT), (3.13) 

K=1 ,3 ;5 , * • • 

vihere and are unknown functions of 5 . 


Substituting y^ from eqte-tion (3.13) into equation (3.4), one 


obtains 


r r , /. \ ,K^2 2 - 




o ir 9 9 ^%r 1 

. . (u^Y.^) ^-(|)V ^> 003(1 KOx) 

as a§ 




“Bu (^Qt) 




+ {-(“) s V 37 + 2 6"o'‘’37’ 

M (f ^ '¥ “-) 
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2*2 *2“ 
“o* 


4 


P {sin ( t) + sin (^~ J2 t) } 

dr 2 


2.2 ^2^ 


u 6 
0 


■K 


4 2 

^ dr 


{cos (SlinO + cos (^nt) }] = 0 


(5.14) 


As stated earlier, the regions of principal primary instability 
can be obtained quite accurately by truncating the series in equation 
(3.13) only at K=1, Thus, equating coefficients of sinC-^S^i) and 
cos (^Qt) from both side of eqi^ition (3.14), one gets 




2 ..2 

u 6 


dC 


i + ( Uq +Y + -^ - u^6) 


1 2- 

^^2 “ 4 ^ ^ ■'^^o ^dS 


= 0, (3.15a) 


and 


d^y 


2 2 


dC 


^ + (u2 +y+-~- + u^6) 




d?' 


112- 
2 - 4 Y^+ 3Uq£5 


= 0. (3.15b) 


It should be noted that equations (3.15a) and (3.15b) are 
coupled differential equations and their solutions can be written as 


'J'i(5) = , i = 1,2, 


(3.16) 


vhere and ^2 " * 

Substituting equation (3.16) into equations (5.15), one gets 


2 2 


{X'^+(U^ + Y+ - u^d) C?) - Su^ BX ^2 = (5.17a) 

uV 

{x'^+Cu^ + y+ -2~- + u^d) x^ -J a } +BU^ nx = o. 


2 


(3.17b) 
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Por noa-trlvial solutions of Q and C , 




-e^Q fix 


^0 


2 2 


( 3 , 18 ) 


Prom equ^-tion (3.18), ei^t values of the X*s are obtained. 
Hence, the general solution of eqtiations (3.15) is given by 

® ¥ 

’"-i = I ® » i=1»2 


(3.19) 


Equations (3.15) hold good for any span of the pipe, She 
solution for ar^r span of the pipe, say the nth span, can be written 


n , , ? n Xj5 i=1,2 

t, (5) = I 0,. e ■* , _ „ 


n=l,2,...,H 


( 3 . 20 ) 


Ihe unknown coefficients C. . can be determined -using the 
boundary conditions imposed by the supports. These boundary con- 
ditions are as follows: 


(i) for zero deflection at all the supports 


( 0 ) = 0 , 


(3.21a) 


( 1 ) 


(3.21b) 


for i = 1,2, 
and n — 1,2,.*., H 
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(ii) tor zero moment at the end supports 


= 0, ■ (3.21c) 

= 0, (3.21d) 

for i = 1,2 

where primes denote differentiation with respect to 5 . 


(iii) for continuity of slopes and moments at each intermediate 
support 1 

= 'f' (O) , (3.219) 

(l) = (o) » (3.21f) 

for i = 1,2, 

and n = 1 ,2,. .. ,M-1 

Moreover, using equ&tion (3.17a), 0.. and C„. can be related as 


G. . 
ID 


0 

1 * 


^.1 


" ^p-i» 3— 
>4 


,.,8 


Q 


( 3 . 22 ) 


Substituting equations (3.20) and (3.22) into eqr&tions (3.2l), 
one obtains 811 homogeneous equations in 8H urknowns. Itor non-trivial 
solutions of the constants consequently °2j ’s), determinant 

of the coefficieris of the resulting homogeneous equations must varhsh, 
Erequencies giving zeros of this deteimnant define the bouMaries of 
liie primary instability regions. 
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ETow, t;o deteriaiae -che regions of secondary insftability, the 
displacement, is expressed as 

y (C,-r) = I iL(5)sin(ij2T) +y (5) cos (Inr) (3.23) 

^ £=0,2,4,... 2r K 2 

Substitution of equation (3.23) into e^tion(5i4) results in 
equation (3.14) vdth summation now being over K = 0,2,4,ii. . The 
regions of principal secondary instability can be predicted with 
reasonable accuracy by truncating the series in equation (3.23) at 
Z = 2. Equating the constant term, and the coefficients of 
sin (fir) and cos (S^t) terms from both sides of the resulting equ^ion, 
one gets 
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Substituting equation (3.25) into equations (3.24), oce gets 


yi Q U 0 ^ 

{X + (u^+ Y+ As-) X } C + ^gu^ saxc + u?6X^ 


0 2 ^ ^o'^^ 0^=0 


.4, , 2 


O N ,2 ? 


{X + (u^+Y + “j— ) X ^ [f } G^- 2Bu^fiX C^-B u^6fiX C^= 0 

2.2 

<4- 2 , *^o 2 2 2 

{X + (u^+Y + —g" )x - Q } c^+ 2 ^u^QX C^+ 2u^<SX G^= 0 


(3.26a) 

(5.26b) 

(3.26c) 


Por non.-tri-vial solutions of C^, C^, and 



j 5S)X 

2 r,2 

u 6X 

0 

~ Bu^ 6GX 


-2Bu ox 
o 

2u^ 6 X 

0 

2^q nx 

X +(Uq+Y+ 


(3.27) 


Irom equation (3.27), twelve values of the X’s are obtained. 

It is obvious that two values of the X's (say the first two valises) 
are zero, . So general solufcion of equations (3.24) can be written as 


’I' =C + g’ 5+ I C„. e^' (3.28a) 

^o o1 02 oj 

. 12 X.S 

i|). == y G . . 6 , i = 1,2 (3.28b) 

0=3 ^ 

As implied from equations (3.26b) and (3.26e) the coefficiecis 


. 12 X.S 

fi “ I “ij ^ 
j=3 ■' 


C ard. G corresponding to zero values ot the X’s are zero. Por 
1 2 
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remaining ten -values of the X.’s, C and C can be expressed 

D iD 23 


in terms of using equations (3.26) . 

How, the solutio'n for any span of the pipe, say the nth span, can be 
written as 


12 


n ^ rP- X rP ^3^ 

*0 = <=01 *° C 2 I 0"^ = 


3=3 


03 


12 - 


,b- v ''j' 

^i = I ® , 1 = 1,2 

3=3 ^ 


(3.29a) 

(3.29b) 


for n = 1,2,.*.,F 


The unknown coefficients C .'s can be deteimined, using the 

boundaay conditions given by eqijiations (3.2l)* Of course, the 

subscript i will now have three -values 0,1, and 2. Substitutin g 

equations (3.29) into the bouMaiy conditions, and equating the 

determinant of the coefficients of the resulting equation to 2»ro 

(fot- non-trivlal solution of C .’s), gets the frequencies defining 

^3 

the boundaries of the secondary instability regions. 


3.6 Third Method 

In Chapter 2, natural frequencies of a periodically supported pipe 
were obtained using the propagation constants of free hazmonie m-ves 
travelling in a similar infinite structure. The same concept can also 
be used to d^erralne the regions of instability discussed in pre-vious 
sectioacF. It has already been seen that frequencies, fi»e,at which 
equations (3. 15) are satisfied with appropriate boundaiy conditions of 
an lif-span pipe, give the boundaries of -the piimaiy instability regions. 



Likewise, frequencies at which equations (3,24) are satisfied with 
appropriate boundaiy conditions of an pan pipe, give the bouBdarLes 
of the secondary instability regiom. 

Comparing equations (3,15) with equation of free vibration of 
a system, frequencies O' s (defining the instability regions) can be 
thou^t of as the natural frequencies of a periodically supported 
pipe, motion of which in each bay is governed by coupled differential 
equations (3*15)'^ • Use of ware ajproach to determine the boundaries 
of the instability regions is presented in this section. 

Let us assume the pipe imder consideration being infinite in 
length, first of all, we discuss the primary instability regions 
only. For K = 1 approximation, using eqi:e.tion (3.13), the displace- 
ment of the pipe in apy bay, y^, can be written as 

y^ = sin (“Ot) + cos (l-nT)'; (3.30) 

vdiere X^ and satisfy equations (3,15), Since these are coupled 

diffearential equations, dispLacemerJt y^, in general, will contain both 
sin (■”• 52 t) and cos terms. We associate two differert propagation 

constants to these two terms. As it will be shown later, for the 
determina-cion of natural frequencies (i.e,, defining the regions of 
instability), these propagation constants iceed to be calculated xmder 
restricted conditions. Considering ary span of the pipe, say AB 
(Figure 2, l), for periodic structures, one can write 

t Hereafter, these frequencies, £2 , defining the boundaries of the 
instability regions will be referred to as natural frequencies. 
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6-r^ + 8-p^^ cos 


^Bx 




®Ax®^^ ^ cos(Jot), 


A7 


(3.31a) 


and Mg — ®An 005(2 Or) 

— Mx 1 -.1% 1 

= M . e sin(— fhr) + M. _e cos (-Ht), 


(3.31b) 


where 6 aod M lepresent the slope and the bending moment (noi>- 

dimensional) with subscripts A and B referring to the eMs A 

and B and subscripts x and y referring to the contribution of 

the sine and the cosine componerts, respectively. ^ and y are 

X y 

the propagation constants a.$sociated with the sine and the cosine 
waves respectively. 

As explained in Chapter 2, to determine the natural frecpeccies 
of an !K-span pipe, aty dynamic distuibaruje is considered in terms of 
two opposite going waves. 


let 


and 


= M sin(~fiT) + M cos(-^£br) , 

^ tr ' ^ 


x"*" 


= M sin(4fi'^ + 1 _ cos(4£^), 
2 ' y '2 ' 


(3.32a) 

(3.32b) 


be the bending moments at the first support due to the positive 
going and the negative going wave's respectively, IQie respective 


bending momeits due to these waves at the last support will be 

Up*' . _ ^ 

M e ^ sin(^fiT) + M e ^ cos(-r£2T) 


(3.33a) 


r 
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M = M e * sinC'rfJT) + M e ^ cos(^t) 

ir yT ^ j- 2 


(5.35b) 


Since the total bending mcments at the end supports naist vanL^ 
(for these are simple supports), one gets 


= 0 

(3.34a) 

= 0 

(3.34b) 


V+ *Sr =° 

Substitutin f equations (3.32) and (3.33) into eqi;&tions (3.3.4) 
and equating the coefficients of the terms sin('|^JT) and cos('^t) 
fi^om both sides, one gets 


M + M =0 


M + M = 0 

y" 


ITP 


Up 


M e ^ + M e ^ 


Up _ UP 
1 e ^ + M e ^ 

y+ y- 


= 0 


(3.35a) 

(3.35b) 

(3,35c) 

(3.35d) 


Equations (3.,35) are satisfied by either of the following 
coniitions 

Up’*’ Up" 

(e ^ — e ^) =0 with M = M =0 , 

y+ y 


UP mi,, 

y » y^ = 


(e ^ - e = 0 with 


= 0 . 


(5.36a) 

(3.36b) 
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Sea. Gupfca f 5 1 3 has shown "that; for the d.etemuLrja'tion. of natural 
frequencies, only propagating waves (with = o) need to be considered. 
So equations (3,36) reduce to 

iUll' iSli? 

(e - e = 0 with M = M =0 (3.37a) 

y" 

iUV^ iNp? 

(e - e = 0 with E = M = o (3,37b) 

x+ X" 

As usual, the values of I p? _ pT 1 and Ip'*’ - u~ I in 
p3ropagation bands will vaiy between 0 and 2 it , Hence, equations 
(3.37) can be simplified to 


and 





d=Oj 


1, 2, with M = M = 0 (3.38a) 
y+ y- 


3=0, 1,2,.,., H with i = M =0 (3.38b) 

X’*' X** 


How the grafhical approach presented in section 2.3 can be ptaed 
to find out the natural frequencies, ELrst the curves of p* (with 
1 4. = M _ = 0 ) versus 0 and p* (with M , = M = O) versus ft 

y y y x+ x~ 

are to be drawn, Ihen, by dividing the range of p* in the propagation 

bands into H equal divisions, one gets H natural frequencies after 
deleting the hipest frequency in each propagation band. Bie two sets 
of frequencies obtained in this manner give the boundaries of the 
primaiy instability regions. 


The propagation constants with 

and p with M = i = 0 (i.e., M = O) 
y x"*” x“ ^ 

discussed in the f ollowing section. 


1 ^ = M =0 (i.e., i = 0) 
y+ y- y 

can be determined as 
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/ 

3.6,1 Detenoimtion of Itopagation Oonstants 


To determine the propagation constants, receptance method is used. 

In case of the propagation constant v with M = 0, equation {3.31b) 

y 

reduces to 

\ sin(|j2T) (3.59) 


It should be noted that due to coupled differential equations 
(3.15), receptance g. , in general, will be of form 

J-J 

vsfaere 8 and . are the sine and the cosine components of the 
slope at end i due to a non-dimensional unit harmonic moment 
sin(^J7T) at end j, 

Ibcpressing slopes in equP^tion (3.31a) in terms of the moments 
and throu^ the receptances, one gets 



1, + 

ibc 




(ef, M. + 

AA Ax 


V V '2"^^ 

^ . e'^ (0^ 5^+ 0^ Mg^=o=(^) 


(5.40) 

Substituting equation (3.39) iit:o equation (3.40) and equating 

-| ^ 

the coefficients of the terms sin(~0T) and cos(-^t) from both aides 


of the resulting equation, one gets 


2u 




- 


(3.41a) 
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gy _ n 

■' ~ ^BA ~ ° 


(3t41b) 


Erom equation (3.41a), two values of (with 1^=0) are 

obtained, one for the positive going wave and the other for the 
negative going wave. Equation (3.41b) gives the value of V vibldx 

y 

is of no use for the- present computation, as we need to calculate the 
values of y with = 0* The receptances to be used in 

j 

equations (3#4l) are determined in the manner discussed below. 

lo deteCTine the receptances and apply a non- 

dimensional unit harmonic moment sin(^ fix) at the end A of the 
span AB (ligure A-I.l). Solutions fcsr eqi;e.tions (3.15) are now given 
by equations (3.19)» with the following boundary conditions: 


X^(o) = 0, 

Y/O) =0 , 


x^(l) = 0, 

O 

II 

(3.42) 

~x!|'(o) = 1, 

-Y''(0) = 0 , 


—X” ( 1 ) *= 0 j 

« 

O 

II 

V 



Use of equation (3.1,9) into equations (3.42) yields the following 
system of eqv&tions 


ZC = 


( 3 . 43 ) 
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where Z= 






-x; 


-x^ > 




.X?t 


:i -1 


2 h 


2 2 


-X, 


-X^a'5 


V 

-sS 


'8 


2 

e^® 
-Xq e 


8 


-8 , 

J 

V 

■^8'^8 


■x|v'^ -XjV^^..- xlv® 


2;i2 

U 0 

o 


"■n + K +1' + — -^o*) 'n- — 


t = 
n 


^ xi=^1 ^ 2^ • • *3 9 


eu^ax^ 


°° ^°11 °12 °13 °14 °15 Se °17 *^18^ ’ 

V t 

T ={0 0 1 0 0 0 0 0 }, 

1 

and t denotes the transpose 

Knowing G^^’s fiom equation (5.45), G^^'s can be deteCTiaed 
using equation (5.22). Hie recejtances 6^ and can be 

expressed as 


and 


(5.44a) 

(5.44b) 
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8 


viiere 

3^ 

M 

= x»(o) = 

I 

n?=1 

°1n 



^AA 


8 




= y«(o) = 

I 

10=1 

°2n 






8 




^BA 

= x»(i) = 

I 

n=1 

°ln 

X e 
n 




8 



and 

BA 

= y;(i) = 

I 

ia?=1 

°2n 



To calculate the receptahces 3^ and apply a noa-diioensiODal 

unit harmonic moment sin('^{^T) at the end B of the span AB» Hie system 
of equations obtained in this case is 


ZC = , (3.45) 

vfaere T» = {0 0 0 1 0 0 0 0}^ , 

and t denotes the transponse. 



(3.46a) 

(3.46b) 



and liie are solutions of equation (3»45) with the obtained 

by using equation (3.22). 

In case of determining u with M = 0, equation (3.31b) 

j X 

reduces to 

Mg = cos(^Jh) = cos(is2T) (3.47) 

Ihe slopes in eqiation (3.31a) are expressed in term of the 
momeobs and throu^ the recepbances (now to be obtained 

with a nort-dimensional unit harmonic moment cos(^nT)). Ihen, 
equation analogous to eq\;ation (3.41a) is obtained as 

els ® - sL = ° 

Itiile calculating the recepbances with a non-dimensional unit 
harmonic moment cos(^nT)^ vector in equation (3.43) and I’ in 
equation (3.45) will be 

= {0 0 0 0 0 1 0 0 }* , 

1 

and T» = {0 0 0 0 0 0 0 1 }^, 

respectively. 

Now, to determine the regions of secondary instability, using 
equation (3.23), the displacement of the" pipe in any bay, y^ , is 
written as 


y^ = Yq + Xg sin(nT) + Y2 cos(Qt) 


(3.49) 
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Proceeding simLlarly as in "the case of primary instahilitiesj 
equations analogous to equations (3.41a) and (5.48) can again be 
obtained for* propagation constants. The unit harmonic moments to be 
used} for calculating the receptasaces, in this case win be sin{DT) 
and cos(f1 t) instead of sin(~!i2'^ and cos(^t) j respectively,. 

3.7 J5ffect of Neglecting Mass Eatio Earaiaeter 


If the terms with mass ratio parameter, & , are neglected in 
equation (3.2), all the methods for deteiminit^ the instability 
regions are simplified considerably, !l3ae modified equations, viien 
the terns with 6 are neglected, can be obtained by letting g = 0 
in the analyses presented. It shorld be rwted, however, that m^ 0 
is not implied by letting g equal to zei*o. 

In the analysis of the first method, if g is set to zero, then 
for K = 1 approximation, equation (3.10) reduces to the foUtwdng , 
form 


Hence, the boundaries of the primary instability regions can be determir^ 
from the eqmtions 

2 2 

2 ^ u 5 _ 

[I] + |-s]+ (u2+y+-|— -u;6) [r3=0, (3.50a) 

aoa Il22=-4 

Eqt;ation (3.12), giving the secondary instability regions for 
K = 2 appr-oximation, reduces to 
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•^1 

^13 

22 

D 

31 

D 

33 

o 

11 

and 



*11 

Di3 


^1 

®33 


= 0 


= 0 , 


give the boundaries of the secondary imtability regions. 

In the analysis of the second method, discussed in section (5. 5 _ )» 
if 3 is ne^ected, then equations (3.15) reduce to the following forms: 


dV u^d^ d^X 2 

+ (u^+ Y+ — — - u^6) i - % 


= 0 


(3.5la) 


l.iL .0 

dt d5 


(3.51b) 


Thxis, we get uncoupled differential equatioce for and Y^. 
Bie solutions for equation (3.51 ) can be written as 


S ‘ J. °i 


^in^ 


n=1 


n 


4 

I 

n=1 


’^l ■= I °211 ® 




X *s and X,^ 's are the roots of the polynonnal 
in 2n 


(5-.52a) 

(3.52b) 
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2 2 


^ 2\2 0 
\ + (uq +Y + — -U^5) X =0 , 


2 2 


X-* +(u^ x^.i!. „.o , 


2 ”0“^ ~ 4 


respectively, Ihe constants and can be determined usir^ 

the boundaiy conditions given by equations (3,21). Substitution of 
equations (3.52) in the boundaiy conditions and for non-trivial solution 
of (and the determinant of the coefficients of the 

resulting equation is equated to zero. In this ease, we get two 
determinants, each of size 4ir, zeros of ifiiieh give the boundaries of 
the primaiy instability regions. 


Itor the secondaiy instabilities, if g is neglected, equations 
(3.24) reduce to 1h.e following foims: 


1 ? 


(uf +Y + 


0 




2 - 

-0X2 


= 0 


(3.53a) 




+ Y + 




2 

+ u„6 

0 


dY 

2 

dC 


= 0 


(3.53b) 


4- 

d Y, 


de 


T * 


2 - ’ 2 

)- 4 - 


= 0 


( 3 . 53 c) 


d 5 


One boundaiy of the instability regions is obtained fiom the 
solution of equation ( 3 . 53 a) and the other from the solution of 
coupled differential equations (3.53h) and (3.53c). Solutions for 
equations (3*53) can be written as 



82 








Y = I c . « 23 


° 3“, -03 ° V^osf- 


aM Y = I 0 
3=1 ^ 


^23 « 


viiere are the roots of the polyrKsnial 


2 2 
u 5^ 


+ (u2 + Y+ ) X^ ~ Q =0 


and Xp.’s are the roots of the polynomial 

2 2 ’^o^^ 2 

x2 + (u^+ Y + —y) u^6 


2u^6 

o 


nV 

4. 2 ^^0° \ 2 2 

X +(ug+ Y+ --y) X - ft 


= 0 


Ooefficdenfcs 0^ ’s, C ’s and 0. ’s can be detenaLned 
In ' on 2n 

iising the boundaiy conditions of an H-span pipe. Eroceediig 
similarly as in section 3. two detenaLnants are obtained in this 
case, !Che zeros of liiese detenninants give the boundaries of the 
secondary insbability regions. 

If 3 is ne^ected, the mve approach for liie determination 
of natural frequencies developed in section 2,3 be used directly. 
With B= 0, one gets unccxipled differential equations for ^ and 
f , One boundary of the instability r^ons is obtained sol-ving 
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for and the other is obtained from the solijtion of . As 
Chapter 2, assume the solutions for aM Y^ of the form: 

X = T 0 


XU 


’l = J. 




a=1 


2n 


vtoere liie X^^^' s and the are the roots of the polynomial 

,4 /:2 




and - (u^+ Y + 

0 


“o«' ^ 2 1 ,2 fi" „ 

- 2 “^ V) >^-- 7=0 


(3,54a) 

(3.54b) 


respectively. Using equation (3.54a) instead of equation (A-1,5), 
propagation constants can be determined followir^ the method outlined 
in section 2,2,1, Erom these propa^tion constants, natural 
frequencies can be obtained using the grajfaical method discussed in 

4 ' 

section 2,3, Ihese frequencies give one of the boundaries of the 
instability region. Simileirly, the other boundaiy can be obtained 
by using equSition (3,54b) instead of equation (3,54a). 

Eor the secondary instability also, one of the boundaries can 
be obtained in the manner described in the case of primary iretability, 
only equation (3,53a) has to be used instead of equation (3,5l), 

!Hie other boundary is obtained from the solution of coupled differential 
equations (3,53b) ard (3.53c), Hie method developed in section 3,6 _ 
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has bo be used "to detennine this boutidaiy. It should be noted, however, 
that in the present case, the order of the characteristic polyrtoiriial 
equation for x' s is reduced to only 6 (from 10 in section 3.6 ). This 

reduces the computational time considerably. 


3.8 Effect of lamping 

If the damping in the pipe is included, equation (3.4) can be 
written as 


4- 2 p ^2- 2- 

(l+a-^) — ^ + {u (l + 6 cosOt) +y}-— ^ + 2ru.(i+6 cosQ 

at O Z O 

^ 3C 


35 


2 - 


- 0 u 6Q sin OT -2^ + q -^ + = 0 

0 35 car ^^2 


(3.55) 


, - - r I -,1/2 E* 


% = 




{(m^+m )EI } 


J72 


E'^ is the loss modulus of the pipe material, and 0^ is the 


viscous damping on the motion of the pipe. 


In the following analysis, hysteretic type damping has been considered, 

I 

i.e^, afi = is replaced by , where q^^ is called the loss factor 
of the pipe material. 

The regions of parametric instabilities of equation (3.55) can be 
found out by determining the periodic solutions of period 2T and T. 

All the three methods discussed earlier can be used to compute these 
regioiis. The dianges necessary in the formulations described earlier 
for undamped case, are given below. Rest of the procedures remain same. 
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If the first method is used, then the teme D „ and D in 

12 21 

equation (3. 10 ), giving the regions of primariy instability, becoiie 


I) 


12 


= \ [s]--1ti n [i]-gu Q [p] , 


aod = I [S] [l] + g u^Q fp] , 


respectively. 

For the regions of secondary instability, inclusion of damping 
changes the terms in equation (3.12) as follows : 


^23 = -hb [ S] - 26 Uq Sfi [P]-n^fi[l], 

^2 " % f ^ ^ ’ 

and the other terms remaining same as defined after* equation (3.12). 

In the second method, inclusion of damping teisas charges equations 
(3.15) (for the primary instability) as follows : 




M 2 2 A d^Y dY 


d'^'Y 


2 *2 
U 0 


d^Y. -2 


Y^ =0 




(3.56a) 


dX. 


1 , (,2 , ^ , ,2 ^ . I ^ . e ^ 

(3.56b) 


4 ' 0 Z o ^ 1 

dC U ^ d5 dt 


+ =0 


Solutions for equations (3.56) can be determined as discussed in 
section 3.5__. Phe characteristic polynomial given by equation (3.18) 
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now beeomes 


2 2 

An U 5 2 

+(Uo + Y + 6) 


1 ,4 1 

- 2n^ X - auQx-gn^jJ] 


1 ■> 4 1 

o + P'^six + -r nn 


,4/2 ^ 2 V 2 £i^ 

X + -Y“ + 5) X'" - — 


(3.57) 


For the regions of the secondaiy ii^tability , equatiocB (3,24) 
reduce to the following forms : 


6^ d^ 


^0^/2 0 " , ^0 2 ^2 1 

^ + (u^ + Y + __) __ + u 5 + 1 g u 5n ^ = 0 , (3.58a) 

; dC dr ^ 


*^"^^2 2 ’^o 2 - ^ ^2 '^^2 

^.(Uo^Y.— )^-. 




(3*58b) 


4-* 2 

,2 


^ d^ 


r -%-) 7 ? - ^2 + "b + 26 U^a ^ + b a 1 


2 ' 2 
^ d5 


^ ‘-dS* 


o dc C- “2 


2 

+ 2 u^ 6 ■;— = 0. 

o aC 


(5.58c) 


Solutions for eqmtions (3.58) can be obtained as discussed 
earlier by changing the chaiacteristic polynomial, given by eqiation 
(3.27), to 



87 


2 2 
^ 

x'^+(Uq+Y+-^) -J 


2 * >,2 
u 6 X 
o 


— gu^finx 


2U=SX 


• 4/2 0 \2>i 4- 

X +(uq+y + — ^^)x -Ji^x -agu^nx - ti^a 


\ X + n^A xV%+ Y+ 


= 0 


( 3 . 59 ) 


The ■fctiird meliiiod discussed in ssciiou 3.6 , can be used as 

vdaen damping is present. The only difference being that the X*s for 
the primary instabilities are to be calciiLated by using equation (3.57) 
and for the secondary instabilities by i;^ing equation (3,59), 


3.9 Results and Discussions 
3.9.1 Computations Iterformed 


Using the analyses presented in the previois sections, the 
following computations have been performed. The regions of the primary 
and the secondary instabilities associated with the first two modes of 
a two span pipe are determined. The instability regions have been 
presented as plots of the non-dimensional frequency, Q versus the 
excitation parameter, 5 . Values of 6 upto 0.5 have been considered in 
all the cases, 

(i) The regions of instability have been obtained with = 2, 

Y = 2 and P = 0,5. All the three methods have been used for the 
purpose of comparison and cross checkirg, 

(ii) The effect of the mass ratio parameter, p , has been studied by 
taking various values of 8 like 0,2, 0.6, and 0,8. The other 
parameters being u^ = 2, and y - 2, 
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(iii) Results have also been, obtained by neglecting g altogether for 
^0 ~ Y = 2, These results are then compared with those 

obtained by taking $ = 0,5, other parameters remaining the same. 

(Iy) To study the effect of the constant njean fluid velocity on the 
regions of instabilities, results have been computed with 
■U-Q = 0»5, 1.0, 2,0, and 2.5. The other parameters being Y = 2, 
and 6 = 0.5. 

(v) The effect of the fluid pressure has been studied with values of 
Y = 0,1,2, and 4 1/^1110 mainbaining other parameters unchanged, 

(vi) To study the effect of damping on the regions of instabilities, 

domputations have been carried out by taking viscous damping 

parameter, n = 0.2, 0,5, aM 0,75. The other parameters being 
c 

Uq = 2, Y = 2, and 3 = 0. The results have also been obtained for 
various values of tbF hyste retie damping, n^« 

The results mentioned in paragraphs (ii) - (v) have been ctotained 
by using the wave approadb outlined in section 3.6 . However, to 

study the effect of damping, numerical results have been obtained using 
all the three methods, 

3.9.2 Comparison of the Methods 

The regions of the primary and the secondary instabilities 
associated with the first two modes of a two span pipe have been computed 
using all the three methods discussed in section 3.4, 3.5, and 3.6, 
respectively. 

Ror the first method, equation (3. IO) has been used to obtain the 
regions of the primary instability. The comparison functiore ^^sed 
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to evaluate the matrices P, Sj R are the normalised mode shapes of a 
two span beam. Regions of the secoodaiy instability have been obtained 
by using equation ( 3 . 12). Six mode approximation was found to be 
sufficient for the instability regions associated with the first two 
modes. 

To get the results by the second method, procedure discussed in 
section 3.5 has been used. 

Ydiile using the third method, the propagation constant v with 

X 

= 0 , is determined using equation (3.41a). Thereafter, natural 

frequencies are obtained by using the grajhical approach. These 

frequencies give one of the boundaries of the instability regions. 

Similarly, the other bounda]y is obtairted by determining the propa^tion 

constant p with M =0. The regions of secondary instabilities have also 
y X 

been obtained as discussed in section 3 . 6 _. Por a two span pipe, 
results for both the modes are obtained frcm the firfft propagation band . 

It should be noted at this stage, that the natural frequencies 
of a two span pipe can be classified into two sets. Values in one set^ 
consisting of the odd-numbered natural frequencies are given by the 
natural frequencies of a single span 3 a. nned- pinned pipe. The olher 
set defining the even-numbered natural frequencies consists of the 
natural freqxKrffiies of a single span clamped-pinned pipe (see Section 
2.5.3 ) . The determimtion of the instability regions ultimately 
reduces to finding the natural frequencies. Hence, all the results 
for the instability regions of a two span^can be c^tained from, those 
of a single span pipe with appropriate end conditions. Hcsjever, this 
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simple extension from the cases of single span pipes is not possible 
for all the modes if the nuniier of spans is more than two. 

Results obtained by all the three methods have been plotted in 
ligures 3* la and 3*1b, The results for the primaiy instability have 
also been presented in Table 3-1. It is seen that the results obtained 
by all the three methods are in excellent agreement v/ith eadb other. 

The relative merits and den^rits of eadi method are disctBsed 

below. 

The first method, thou^ simple to use, has certain drawbacks. 

Ror example, on evalmting the matrix [ P ] , one finds that the diagonal 
elements p^^, Pggj etc,, are zero. Thus, for one mode approximation, 
the elements ^21 equation (3.10), and the elements D^g, 

Rgj, and iu. equation (3.12) become zero. This readers equation 
(3.10) and (3.12) independent of the value of g , Thus, the effect of 
the Coriolis acceleration in the equation of motion (3,4) is not accounted 
for. Moreover, it is obvious that parametric instabilities are obtained 
ohly in the modes used in approximating the pipe displacement. To use 
this method, the mode shapes of an H^-span beam must be first obtained. 

Hence, it becomes tuiwieldy especially for large nusher of spaans in the 
pipe. 

-fine 

As compared to the first method, ^second method requires oo preidous 
knowledge of the mode shapes of an IJ-span beam* Moreo^rer, the iiBtability 
regions associated with all the modes are obtaioed simultaneously by 
determining the zeios of the same deteiminant* However^ 14ie order a£ this 
determinant increases with the nunfoer of spans in the pipe* 
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Values of frequencies defining the boundaries of the primary Instability regions associated with 

the first two modes of a two span pipe# 
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As compared, to the other methods, the third method is most suitable 
for large number of spans in the pipe» Once the propagation constants 
are obtained, the instability regions for any number of spans can be 
obtained by dividing the ranges of \i* ani y* (in the propagation banis) 
into the same number of eqi:a'l divisions as the ruxnber of spans in the 
pipe. Thus, in the third method, the amount of computation is independent 
of the number of spans in the pipe. 

3.9.3 Effect of the Mass Batio Ihrameter 

Tigures 3- 2a and 3.2b shew the regions of the primary instability 
associated with the first two modes of a two span pipe for various values 
of the mass ratio parameter, P , It is seen that increase in 3 sli^tly 
shifts the instability regions to lower frequencies. However, the ^ift 
reduces with increasing 6 . The maximum shift in the regions is 
approximately within 3% over the range of g given by 0.2 ^ B 5. O*®* 

Hgures 3,3a and 3.3b show the regions of the secondary instability 
associated with the first two modes for varioie values of B . !Qae effect 
of variation of B on these regions is same as that on the regions of the 
primary instability. 

As stated in adapter 2, the negligible effect of B ^ould mt, 
however, be interpreted as negligible effect of the mass of the fluid, 
This is because the non-dimensional frequenqy n also contains the 
term m^. Thus, the narginal effect of the value of B only signifies that 
the terms containirg B » liEe the Coriolis term, have little effect on 
the plots of 52 versus 5 . 
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^•9*4 Effect of Ee^ecijing the Ifess Batio ihrameter 

To study the effects tie^ecting g altogethexf results have been 
obtained using the third method discussed in section 3.7. Eigure 3,4 
shows the regions of instability associated with the first mode of a 
two span pipe with u^ = 2, Y = 2, arid g = 0, Eor the purpose of comparison, 
the same regions obtained with u^ = 2, y = 2, ani e = 0.5 are also shown 
in the same figure, it is seen that the results obtained by neglecting g 
altogether have an accuracy of the order of 9^ for the priisaiy 
instability region. The legions of the secondary instability are not 
disting-uishable for g = 0 and g = 0,5. However, as slated already in 
section 3.7, the computational effort reduces considerably if H is neglected, 

3.9.5 Effect of the Constanb Mean Velocity 


The method discussed in section 3.6 has been used to study 

the effect of the constant mean velocity, u^, on the regions of instabilities. 


Figure 3.5a shows the effect of u^ on the legions of primary 
instability associated with the first mode. It is seen that the regions 
become wider and are lifted to lower frequencies with increasing u^. For 
large value of u^, the instability region may even start with 2 ero frequency 
if the excitation parameter is sufficiently higa, Hiis impiies that with tdrese 


combinations of u^ and 6 
coefficient of the term 


, the equi-Talent axial load (represented by Idie 
d^Yi 

— i in equation (3.15b)) is more than the ciitii 
dC2 


load required for buckling. 



Figure 3.5b shows the effect of u^ on Ida® regions of primary 
instability associated with the second mode. The effect is seen to be 
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similar "fco "that Lu tihe first mode. However, the regiotK do not beccsiK 
as wide as in the first modej wilii inciease in u. • 

ligures 3«6a and 3 •6b show the effect of on the regions of the 
secondary instabilities in the first two modes. In these cases also, 

"the regions become wider and get shifted to lower frequencies wilii 
increasing u^. It ^ould be noted that the upper boundaries of the 
secondary instability regions not change mch with increasing 6 , 

Comparing the regions of primary and the ^condaiy instabilities, 
it is seen that the upper boundary of the primary instability shifts to 
hi^er frequency with increasing 6 , vfcereas the upper boundary of the 
secondary instability regions shifts to lower freqisency with increasing 6 . 

As seen in section 3.9.4, the regions of instability do not diange 
significantly by letting g = 0. Thus, the changes in the instability 
regions described abo-ye can also be studied qualitatively from equations 

( 3 . 51 ) and ( 3 . 53 ). The coefficients of the terms having secoM order 

2 - 2 - 

dV dV 

derivatives (like — — , — ^ , etc.) increase with increasing u . The 

dC dr 

coefficients of these terms repacesenb the equivalent axial load on the 
system. Therefore, the ratural frequencies of the systems governed iy 
equations (3.5 I) and (3*53) shift to lower values. 

In equations (3.5l) vtoidi governs the regions of primary instabilities 
the coefficient of the term 1 reduces with increasing 5 , TIhjs the 

natiit?al f!re^encieS| one of the bourwiaiy of the xustahilxly region 

increases to hi^er value with incieasing 6 » Cto. tl^ other handp the 

coefficient of the term increases with increasing 5 . HetMS, the 

dC^ 
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natural frequencies reduce to lower values with increasing 5 , The net 
outcome of these two opposite changes is the considerable increase in the 
width of the instability regions with increasing 6 , 

If we consider equation (3»53) govenaing the regions of the seccsidaiy 
instabilities, the coefficients of the terms with second order derivatives 
increase with increasing 6 , !Ilaus, both the boundaries of the instability 
regions shift to lower values with increasing 5 . As a result, the 
width of these regions does not change as much as in the case of the primary 
regions. Moreover, with S as a fraction, the terms signifying the 
equivalent axial load in equation (3.5 1 ) change more drastically with 5 
than those in equation (3.53). This explains the considerable change 
in the primary instability regions with the increasing value of 6 , 

3.9.6 Effect of the Eluid Pressure 

ligures 3.7a and 3.7b shew the regLors of the primary instabilities 

associated with the first two modes of a two span pipe for various 
values of the fluid pressure, y . It is seen that the instability 

regions associated with both the modes shift to lower frequencies 

with increasing Y . 

ELgure 3,8a aid 3.8b ^ow the regions of secondary instabilities 
for various values of y , The effect of variation of Y on these re gLons 
is similar to that on the regions of the prluBiy instabilities. 

3.9.7 Effect of Damping 

!Qie effect of damping on the regions of instabilities, has been 
studied using the methods discussed in section 3.8, Botii viscous and 
hysteretic types of damping have b^n considered. 
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EigupB 3*93- shows "thG regions of in.sfshili'tics associs-tsd wilSi 

the first mode of a two span pipe for three values of the viscous 
damping parameter It is seen that the presence of damping diminishes 

the extent of the instability regions. A finite critical value of 5 is 
required to have the instabilities. If 6 is less than the ciitical 
value then the pipe is not subjected to parametric instabilities. Kie 
critical value of 5 increases with increasing In other words, 

large excitations are required to excite the instabilities if the dampii^ 
is increased^ The effect of on the regions of the instabilities 
reduces with increasing 6 • 

It is also seen frcm liguie that the damping has more pEroimoaced 
effect on the regions of the secondary instability as compared to the 
primary instability. Par a given value of ’tit® critical value of 5 
is much more as compared to that for the primaiy instability. 

iigure 3# 9b shows the effect of damping on the regimis of 

instabilities associated with the second mode. As compared to its 

effect on the regions associated with the fir^ mode, the damping is 

more effective and the regions become considerably narrower. Por valt^s 

of n jtist more than 0.2, the secondary instability may be totally 

o 

eliminated for 5 even upto 0,5. 

Elgure 3.10a and 3.10b show the effect of the hysfeeretic damidJi^g, 
n^jf on the regions of instabilities associated the first two asides. 

As with viscoTos damping, three values aie a3^ used for 3ii order 

to oompare the effectiveness of the viscous and the l^stesretic "type of 
damping, equivalent values for and are us^ in the compatati«m. 
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This egiii-valence is based oa the same decay rate in the first inode of a 
simpLy supported pipe vihen both the fluid pressure and velocity are aero* 
The effect of ri^ is seen to be more pronouncoi than the corresponding 
value of n^. The secondary instability region associated with the 
first mode ls eliminated totally for =0.191. 5Phe secondary 
instability region associated with the second mode is eliminated even 
for = 0,064. 




CHAPTER 4 


iHV.'iSP IGrAT 101 OP THE EPPECT OP HTHAHIC ABSORBERS 01 THE REG- 10 IS OP THE 
BlR,\I>/IETRIG USTABIIITIES OP A PIPE 

4.1 Introduction 

A dynamic absorber normally consists of a lumped mass that is 

attached by a damped resilient element to a vibrating system, the motion 

of which is found to be excessive, Piequently, the vibrating system of 
fo u nd to 

mass, Mg, is^resonat4_ir' on resilient supports as diown in ELgure 4.1. 

The dynamic absorber is tuned to a frequerKjy close to iiiat of the 
vibrating system. The motion of is then reduced because, in effect, 
the absorber mass, m , is greatly magnified in the nei^bouihood of this 
frequency and controls the motion of 

to 

The djrnamic absorbers have been used successfully to reduce 
undesired resonant vibrations of a systm caused by ’’ground vibrations” 
as shown in Pigure 4.1a. I^namic absorbers may also be emplc^d to reduce 
the force transmitted by a machine to the ideally rigid foundation as 
diowii in Pigure 4,1b. Biese absorbers also find applications in reducing 
the resonant motion of rods vibrating in their longitudinal 3i»des or 
beams vibrating in their- transverse modes [54]. 

In the present chapter, a viscously damped dynamic absorber has 
been tried for a different use from its conventional role of reducing the 
resonant motion. The aim of the present study is to determine whether 
the dynamic absorber can control effectively the regions of parametric 
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instabilities of a pipe. As stated in CJhapfcer 3» in most of the practical 
cases, the pri'ncipal region of instability is of primaiy importance, !Ehis 
is why, in the present work the effect of dynamic absorber on these 
regioas only has been investigated. It has also been shown in Chapter 3 
that the mass ratio parameter, g , has very insignificant effect on 
the regions of primary instability. Hence, in the present study this 
parameter has been defected altogether. Moreover, damping in the pLpe 
has also been omitted. 

In this chapter, first we consider a viscously damped dynamic 
absorber attached at the mid-point of a single span pipe, iliree different 
end conditions for the pipe, viz., pinned-pinned, clamped-clamped, and 
clamped-pinned , are considered separately. The effects of different 
absorber parameters, like the mass, the damping, and the spring constant, 
on the regions of instability have been studied. R)r computation of 
the instability regions, the second method discussed in seeticn 3,5 has 
been used. This is followed by the analysis of a periodically supported 
pip6 spring— mass ijypG absoibsp a*tlJSLGh6d ai IjIis CGJit3?6 of saob span.# 

Til© * W 8 .V 0 apppoacli^ dOYSlopsd. in. sociiion. 3*7 IIS'S bssn iis©d# 

4* 2 Theoretical Ttormulation 

Consider a pipe of length £ restii^ on t3?ansryersely rigid supports* 
A viscously damped absoifcer of mss, spring constant, and 
viscous damping, C^, is attached at the mid-point of the pipe {^gure 4.2). 
Let the velocity of the fluid flowing thrcxigh the pipe, be 

u = Uq (1 + « cos fir), 
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v\3i63!?s is liii© non— dime Qsiooal constant mean velocity j 

6 is the excitation parameter, 

fj is the non-dimensional frequency, 

and T is the non-dimensional time. 

It was seen in Chapter 3 that the effect of neglecting mass ratio 
parameter, 6 , on the princtipal regions of instability was 
negligible, and the analysis becomes much simplified. So for the prei^ot 
study also, the effect of g has not been included. For the regions of 
primary instability corresponding to E = 1 approximtion, the displacement 
of the pipe, y^, is given by equation (3.30) 

y^(C,T) = X^(S) sin (^nr) +f^{c) cos 


where and satisfy equations (3.5 1) 


d'^X 2 ’^o 2 ^ - 

^ + (u^ + Y + — u^ 6) — TT- - — X. = 0, 




4 1 


d^'f 




dC 


1 / 2 0 ' . 2 
^ + (u^ + Y+ — tr" + % 


d? 

I? ■ 


^ — Y. = 0. 


The 


solutions for equations (3.5 1) are given by equations (3.52) 


X.(?) = I 3.^6 
' n=1 ' 

4 ■ ^n^ 

Y^(5) = I ®2n ® » 

' tt=1 


where the X_^»s and the X^^’s are the roots of the podynomials 


in 


X<. (u2 .Y^^-u2 0 x"-r = 0 , 
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aod x'^ + (u^ +y+~ + uj 5) x^ - — . 0, 


respeotlyely. Ihe c,^'s and the o^^'e are unknem coefficients to he 
determined from the boundary conditions. 

4.2.1 Equation of Motion for the Absorber-lfese 


Assume y^ be the displacement of the mass, m . !Ehe equation of 


motion for the mass m is 

a 










(4,l) 


^ = 2 


Substituting the value of y_^ from equation {3.30) into equation 
( 4 . 1 ) and noting nT= yt, one gets 


e.2 


dy Q £0 

"a •*■ °a ar ’'a 5"a = C X,(i) - -f- Y,(i)) t l(k^ 7^(1) 


dt 


C w ±ii^) 

+ -f- e ^ ] , 


(4.2) 


where denotes the imaginary part of the bracketed ^antily. 
The steady state solution of equation (4.2) is 




^ -• 2 k ^ 2k 

m [ 5-- 2 e ] 


lU 0) C ca 

/ ^ a \ . a 

- 4kJ ^ ac^ 


(4.3) 


Let us introduce the following notations 


0 ) = natural frequency of the mdamped absoitier 

a 
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Op — critical damping of the absorber 

= 2m CO , 
a a’ 

c = damping ratio of the absorber 

cx 

C 

a 

“ G ’ 
c 

a = tuning ratio of the absorber 
0) 

a 

^ 0) > 
o 

vdaere is some reference frequency. In the present study, ai^ has 

been taken as the first natural frequency of a pinned— pinned pipe with 

fluid velocity u^ and fluid pressure y . 

Cm 

Using these notations, ^ — can be written as 


0 0 ) 0 

= 2(^) 2(jf) 

(0. 0^ cOp Op 


C CO 


Q n 


( 4 . 4 ) 


XiM ~ 

£l 

Similarly, ^ can be written as 

a 

2 

^a ^ m \ 2 _ / ^Ov 2 /m \ 2 

k “ ^0) ^ ” ^ 0 ) ^ 

Q a n 


~ 2 ’ 
0 
a 


(4.5: 


where is the non-dimensional reference frecpiency 
o 

,“f ^ “p, 2 .2 

= “o* 
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Substituting eijuatlons (4.+) and (4.5) into eolation (4^), 


one 


gets 


y = £I r 
^ ql m l 


) -|y^(|-)) 4 1 (Y^( 1 ) i(l a t) 


4a 0 


1 /'iL^2. 

2 
3 

a 

On simplification equation ( 4 . 6 ) reduces to 


3 (4.6) 


Ja = a Ua X^(-l) + Ba Y^(|)} sin ) + £{a^ ?^{I).BaX^(|-)} cos (-iftt), 

( 4 . 7 ) 


T/diere = 

0 


and B = 
0 


40 0 

Si 

4a 0 

Si 

[1-^(-^)^] (l) 

4a 0 
a 

' " ' " " " '■ " "'f r -r-r r ^ 

r, 1 /^x2-,2' yAx2 

a 


4 . 2.2 Boundary Conditions 

In the present study, the pipes with following end supports 
have been considered. 

(a) 2 Lnned -Binned Bipe 

When both the ends of the pipe are pimed (Slgure 4 . 2 a.), because 
of symmetiy, one needs to consider only half of the pipe and the boui^aiy 
conditions are 

(i) for zero deflection at the left support, 
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y^(o,r) = 0, 

i.e., x^(o) = 0 , 

and Y^(0) =0 

(ii) for z6ro momBnlj af ijli© l6ft support ^ 


(4.8a) 

(4.8b) 


y!,'(o,T) =0, 
i.e.,X!j'(0) = 0, 
and fl|(0) = 0 , 


(4.8c) 

(4.8d) 


vdiGi*© priioes denote differentiation with respect to 5 . 

(iii) IXie to syametry, slope at the mid-point of the span mist be zero 

•') = 0 , 


i-e.r X»(|) = 0 , 


and Ilj(|) = 0 


(4.8e) 

(4.8f) 


(iv) At the mid-point of the span, change in shear force mst be equal 


to the inertia force of the absorber mass, m . Thus, 

'a ^ 


(%) 




^dt^ 


(4.9) 


Substituting y from equation (4.7) into equation (4.9), ore 

obtains 

Sin (IfiT) + fy(j) oos . -1 i)2 [fe^ 5^(1) + 

+ f/^)) sin (i OT) + } ooe (i 0^] , 

(4.10) 
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*ere is known as the mass ratio of the absoiber aal defined as 

m 

r *“ 3» 

a (nu + m )£, * 

P 

Equating the coefficients of the terms sin (|- £h:) and cos (-^ Jhr) 
from both sides of equation (4. 10) , we get 

""e “'a “^^^0 +Bo^l(i» =0 ’ U.Ha) 

( 2 ) + §• fa ^0 5,(4) > = 0. (4. 11h) 

Substituting equations (3.52) into equations (4,8) and (4.II), 
one obtains a system of homogeneous equations 

FP ^ 

vdiere is a matrix of size 8^8 and its elonents are given in Appendix 4. 
The subscript pp refers to the jiQned-pLnrffid pipe, C is a column matrix 

^12 *^13 *^14 ^21 *^22 *^23 ‘^24^^’ transpose. Por 

non-trivial so3ution of C in the above equation, dete rm i nant of P must 

be zero, Erequencies at vdiich the determinant vanishes, give the 

boundaries of the instability regions. It should be noted that even after 
neglecting the mass ratio parameter, 6 , thou^ the governing equations of 
motions (given by eqiiations (3.5 1)) are uncoupled, the boundaiy conditions 
(given by equations (4.8) and (4.II)) are coupled throu^ the damping of 
the absorber. 

(b) Clamped-Olamped Hpe 

When both ends of the pipe are clamped (Figure 4.2S>), then 
a gain because of symnetiy, one needs to ccmider <fflly one half of the pipe. 
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[Che bomdary conditions are same as for the pinned-pLnned pipe (given 
by equations (4.8) and (4,11)) esxcept for equations (4.8c) and (4,8d), 
Since for the present case, slopes at the supports must be zero, 
equations (4.8c) and (4.8d) are modified as 

ylj(0,T) = 0 , 
i.e. , X!j(o) = 0 , 

and Y>(o) =0. 


Proceeding similarly as in the case of a piinned-pinned pipe, 

instability regions are detemined by the zeros of an eighlii order 

determinant |d |, Ihe elements of the matrix D are given in Appendix 4. 
cc cc 

•ihfi 

The subscript cc refersto^clamped-clamped pipe, 

(c) Clamped-Pinned pipe 

ligure 4.2c shows a clamped- pinned pipe with a viscously damped 
dynamic absorber attached at the mid-point of the span. Using equation 
(3.30), the displacement of the pipe can be written as 


sin (^ {h:) + cos (l-fiT), 0 <_ C, (4.12a) 

y^= ^^(Sp sin (^ fit) + 1 ^( 52 ) cos 0 1^2-2 (4.12b) 


\iiere 




I > 

n-1 ' 


(4.13a) 


1 


^ 1 

I <= 


,^2n^1 


(4.13b) 


tt =1 


2n 
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0 ^ 

in ’ 


(4. 13 c) 




4 

I 

n=i 


,^2ia^2 


2n 


(4. 13<i) 


and 5 2 tho nonp-diiiiensioiQal co-^ordlna’fces along thB [tength of the pip^ 

as diown in KLgure 4.2c. The constants o] 's 's, >s, and G^ 's 

In j 2n ’ In ’ " 2a 

can be determined from the following bomdaiy conditions : 

(i) for zero deflection at the supports. S and G 

y](0,t) = 0 , 

= 0 , 


i.e. , xj(0) =0 , (4.14a) 

Y^(0)=0, (4.14b) 

X^(^)=0, (4.14c) 

and Y^(|-) =0 (4.14d) 


(ii) for zero slope at the support E 
y'''(0,T) = 0 , 

i.e., Xlj^Xo) = 0 , (4.15a) 

and ?l|^(o) = 0 . 


(4.15b) 
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(iii) for zero monEnfc at the support G 

»'f) = 0 , 

i.e., X'^' (■^) =0 , (4,16a) 

and Yf(|)=0. (4.16 b) 

(iv) for continuity of displacement at the midr-point of the pipe, F 

i.e.,x](^) = X^(0), (4.17a) 

and y!|(-1)=Y^(0). (4.17b) 

(v) for contimity of slope at the point F 

y!]^(|' > t) = y*^^(0,T) , 

i.e.,Xl,^(-^) = X^^(O) , (4.18a) 

and ?'"'(“)= YJ^(O). (4.t8b) 

(•vi) for contimity of moment at the point F 

ylj’'(^ ,'f) = 

i.e., X!|^(“) = X!,'^(0) , (4.19a) 

and Y!j'^(-i) = Y'j^(o) 


(4.19b) 
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(■vii) ag'din at the point P, change in shear force mrst be egnai to the 


inertia force of the absorber mass, m . !Qiis results in 

8 * 

(~) { + yl"^(o) }= - m — i 

On simplification, equation (4.20) gives 

(A„S](1) H.B„7]ci)} = o 


( 4 ^ 0 ) 


(4.21a) 




(4.21b) 


Using equations (4.13) in equations (4.14) - (4.19) and (4.2l), one gets 
a system of equations given by 


op ' 


( 4 . 22 ) 


where D is a matrix of order 16 x 16, and its elements are given in 

op 

Appendix 4} the subscript cp refers to the clamped- pinned pipe, 0 is 

column matrix with elements 0^2 *^^5 ^^4. ^21 *^22 *^23 *^24 ^1 ^2 ^13 

0^ 0^ (£„ c£, where t denotes the transpose. It diculd be 

14 21 22 23 24 

noted that the order of the matrix D^p is twice the order of D^p and B^, 
This is obviously due to the loss of synmetry of the structure about its 


mid- point. 


For norv-trivial solutions of 0»s in equation (4.22), we get 


Id I - ■ 0 

' op 


(4.23) 


Thus, the zeros of the determinant of the mtrix D^p gi-^ 
bomdaries of the instability regions. 







'\*i (i:.\ '. *v 


,'t ,'N 

>.-v:%-:'. 
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4.3 Undamped I^nainic Absorber Attached at the M3d— Boinfc of Ihch Bay 
of a Periodically Supported Pinite Hpe 


Bigire 4.3a shows a periodically supported U-span pipe. Identical 
undamped dynamic absoibers are attached at the mid— point of each spm. 

!Sie effect of attaching these absorbei« on the regLons of instahilities 
of the pipe is investigated in this section. 


If the absorber is undamped, i.e. , A = 0, equaticn (4.7) describing 
the motion of the absoiber mass, m , in any bay, reduces to 




(4,24) 


riiere 




a 


Let us ccnsider ary span of the pipe, se^ AB shewn in iigure 4.3b. 

Ihe displacement of the piip, y^, in this span is given by equations (4. 12). 

At the mid-point of the span, the change in shear force must be e(jial to 

the inertia force of the absorber mass, m . This result in 

’ a 


(^) [ (^) +y'’’^(o)] 



(4.25) 


Substituting y from equation (4.24) into equation (4.25), we ^t 

£L 


- xy’^(o) ^ Aq x^^(“-)> sin (Iftt) 

+ - Y»''^0 ) + { A^ y;|( 1)} cos {l ax) = 0 


(4.26) 
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Siace this ecjiatic®. holds good for all values of t, one gets 

I 

J x](i) = 0 (A.27a) 

. Yy'2(0) + -J- !? = 0 (4.271)) 

If we examine the boundary conditions for the pipes dealt in 
section 4.2.2, it is found that the boundary ocnditions relating tte 
change in shear force at point F, are coupled in and However, 

for the present ease wilh an undamped absorber, we get uncoupled 
eqi^tlons ( 4 , 27 ) in and lEhus, (one boundary of the instability 

regions can be obtained by solving equation (3.51a) with the corresponding 
bomdary ecnditions. Similarly, the other boundary of the instability 

regions is obtained by solving equation (3. 5 1b) with the correspoadlng 

boundary conditions, 

lo find the natural frequencies of the system governed by 
equation (3.5la), the 'wave approach' outlined in section 2,3 can be 
used, !l!he propagation constants are obtained by using equation (2,3). 

The receptances to be used in this equation are obtained in the manner 
discussed below. 

lo calculate the receptances 3^ and 3^^, apply a non-dim ensioml 
unit hamonic mcment sin (^ £5 t) at the aid A of the span AB, The 


boundary conditions are ; 
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= 0 , 

= ^(0) , 


- = 1 , 

- (^) = 0 , 

= if(o) , 


(4.28) 


(0) , x;"'(l)-x;"^o). 1 r_^ 


4 a “ ^©^1^2 

e of equations (4.13a) and (4.13c) in equations (4.28), results 
the following ^ston of equations. 


in 


% <= = 


(4.29) 


«®e is a matzix (of ord«. 8 x 8) and its alan^ts are glvm la 

Appoadiz 4. 


= g] g], c"' cP -it 

11 12 13 ^4 ^^2 *^3 ^ 4 ^ t 


a?^ = {0 1 0 0 0 0 0 o>^ , 

with t denoting the timspose. 

Ihe receptances 8,, and e_. are given as 


AA 


BA 


v-Vto)= i o' 


£t=1 

4 


in in 


3 r 2 1r/^ 

V - <2’ = 0,^ e 1“ , 


<*ere the 0^^'s and the (^^'s are solitiotis of equatlcn (4.29). 


AB "^BB^ 3- 


lo calculate the reoeptances 3 and 3 

AB j 

■iimaielohal unit hamoole momeat eln (1 8r) at the sad B of the apaa AB. 
The systaa of equations obtained in this case is 
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Z 0=5)’, 
a r 


( 4 . 30 ) 


Til ere I* 


= {000 1 


0 0 0 } 


Hie receptances 3._ and 0__ can be wilttai as 

A3 


1 ^ 1 

= X!,'(o) = I c] X, , 

AB 1 In in ’ 


4 

I 

n=1 

4 


_ o 1 ^9 ^Ir/^ 

Sb ^ ^ ^ ^in ^ln ® 

a=1 


12 

where the Qj^^’s and 0^^’ s are solutions of equation ( 4 . 30 ). 

!Hie other boundaiy of the instability regicais is obtained by 

solving equation (3.5 1b) in the similiar manner as discussed above 

for the solution of equation (3.5 1a). In this case, the receptances 
are to be calculated with a non— dimensional harmonic moment cos (‘ 2 ^)* 


Once the propagation constants have been determined, natural 
frequeacies can be found out by using the graphical approach dev^oped 
in section 2.3. 

4*4 Eesults and liscussions 

4.4. 1 amputations Performed 


Using the analyses presented in previous sections, nmerical 
results, listed below, have been obtained. As in Chapter 3, the instability 
regions have been presented as plots of frequaicy versus the excitation 
parameter. Ihe value of the excitation parameter upto 0.5 has been 


considered 
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(i) The instability regions of a pioned-pinned pipe with and without 
the absorber have been obtained, 

(ii) She effects of Tarious paimeters of the absoiber have been studied 
by obtaining ti® instability regions of a pinned-pianed pipe with 
various ccmbinations of these parameters. 

(iii) Reailts have been obtained to study the effect of u and y on the 

o 

instability regions of a pimed-pinned pipe when an absorber 
is attached to it. 

(iv) The instability regicns of clamped-clamped, and clamped-pinned 
pipes vfith absoibers attached at their mid-points have also been 
computed, 

(v) The analysis presented in secticai 4.3 has been checked by 
determining the instability regions of a two span pipe with 
identical undamped absoibeis at the mid— point of each b^. 

Ctoly a limited range of frequency covering the principal instability 
region without the absorber has been investigated in all the cases. 

4.4.2 Effect of an Absorber on the Instability legions 

To study the effect of a dynamic absorber on the instability 

regions, results have been obtained for a pinned-pinaed pipe with 

and without the absoiber. The method presented in section 3.7 has been 
used for the pipe withcut an absoiber. Using the method outlined in 
section 4, 2, 2^ instability regions of a pipe with an absorber have been 
obtained, Eor the purpose of comparison, both the results have been 


shovn in ligire 4.4 








c’;;- WV-:' 

V-- 


9& 


>. -J'/iV «: 


5> « 'J.. 
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It is seen from iigare 4.4 that the attachment of an absorber 
has maiJced effect on the regicaas of instability. Without the absorber, 
only one region of instability is obtained, vshiie with the absorber there 
are two regions of instabilities over the same range of frequoacies. 
However, the widths of these two regions are quite small as compared 
to the original one. It is also seen that the instability regioas 
obtained with the absorber are distinctly separated from each other. 
Moreover, approximately the entire instability regicja without the 
absorber has been converted to a regiaa of stability when the absorber 
is attached to the pipe, 

4,4,5 Effect of the Absorber Ihrameters 

Effects of various parameters of the absorber on the regicsas 
of instability of a pinned— pinned pipe have been shown in Eigures 4.5 — 

4.8. 

iigure 4.5 shows the instability regions for various values of 
the tuning ratio, a , of the absorber. It is seen that with low tuning 
ratio, the first region is quite narrow and the critical value of 6 , 
required for this region to staid:, is very large. On the otter haad,the 
second region with this tuning ratio is quite wide and critical value of 
6 is rather small. With increase in the first region shifts to 
hi^er frequencies and beccsnes wider. The critical value of the 6 is 
also reduced. Whereas, the second region becoaes narrower and the 
critical value of 6 increases considerably. It is also seen that the 
second region is totally eliminated upto 5 = 0,5 for a > 1,0, Thus, 
the effect of the tuning ratio on the two regions (obtained with the 




absorber) is of opposite nature, 

iigire 4.6 shows the effect of the damping ratio, on the regions 
of instability. The effect is seen to be quite sigaificant. The second 
region of the instability is much affected by the damping of the absorber. 
It is totally eliminated for damping ratios hi^er than 0.05. 

The effect of the damping ratio, on the first region of 

instability can be summarised as follows. Upto a certain value, the 

damping diminidies the extent of the instability regioas, but b^ond 

that critical value, the regions become qrite wide and the critical 

value of 6 starts reducing. This behaviour of damping shows that 

there is an optimum value of ? . This optinum value of ? can be 

a a 

determined only for givoa objectives liie, critical value of 6 i^ould 
be maximum or the width of the regicai should be minimum, etc, 

ligure 4.7 shows the instability regions for a pinned-pinned 

pipe wilh different mass ratios when an mdamped dynamic absorber is 

attached to it. It is seen that with the decrease in the mass ratio, r , 

a’ 

both the regions of instability get closer to each other. Eventually, 
with zero mass ratio, i.e, , without the absorber, only one regicai of 
instability is obtained. Figure 4.8 shows the effect of the mass ratio, 
r , when damping in the absorber is also considered. It is seen that 

a. 

the effect is similar to that of an undamped absorber. But with a 
damped absorber, a critical value of 6 is required to cause instability. 
It is also seen from ligure 4.8 that for low values of the mass ratio 
of the absorber, I , only one region is obtained. For large values of 6 , 
the width of this region is comparable with the region obtained without 
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the absorber (Hgure 4.4), Por large -values of r , two instability 

regions are obtained and the widths of these regions are small in 

ocmparisaa to the regim obtained withcut the absorber (Figure 4.4). 

Shis beha-viour of r indicates that thei«e ezists a critical -value of 

a 

to ha-y-e two regions of instability, JUrtheimore, exaaining 

Hgires 4*7 4.8, it is ccncluded that the damping has more 

prcaiounced effect on the second region (in case it is there) than the 
first region of instability. 

4.4.4 Effect of the Pressure and the Constant Mean Velocity of the Fluid 

Sigure 4,9 shows the effect of the fluid pressure, Y , on the 
regions of instability of a pinned-pinned pipe with an absorber (c = 0.9, 
= 0.1 and = 0.2) attached to it. It is seaa that with increase in 
Y y the instability region becomes wider and the critical -value of 5 
reduces, For sufficiently large value of Y , the3:e may be even two 
regions of instability. Examining lignres 4.4 and 4.9 together, it is 
seen that the second region (in case it is there) experiences similar 
effect of Y as \ ^the first regicxi discussed above, 

Sigire 4,10 shows the effect of the constant mean velocity, u^, 

on the regions of instability of a pinned-pinned pipe with an absorber 

(c rs 0.9, ? = 0,05, r = 0.2) attached to it. It is seen that increase 

in u reduces the ciitioal value of 5 • lEhere may be two regions of 
o 

instability for large values of u^* It is also seen ffom ELgure 4*^ 

that there is no unstable region f or less than If the damping 

in the absorber is incieased to z, =0.1 frcm 0»05j the instability 

a 

region associated. vAth u^ = 1 has been found to be eliminated upto 
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6 - 0,5. Pie in 0 tability regioa for = 2 has already been shorn in 
Eigare 4.4. 

4.4.5 Upes with Other Sid Ctoadlticws 

Besults wilii other aid conditicais indicated yexy siioLlar trends 
of effects of yaricws paraaeteis as observed in the pinned-pimed c£ise. 
Hence, only representative curves are discussed below, one each for a 
(SLamped-pinned pipe and a clamped-clamped pipe, 

iigare 4.11 ^ows the regions of instability of a clemped-pinned 
pipe with the absorber attached at its mid-point. It is sem that with 
an •undamped absorber, two regions of instability are obtained, iHiese 
regions are greatly affected by the damping present in the absorber, 

The effect is similar to that on the instability regions of pinned-pimed 
pipes, She first region of instability is totally eliminated upto 6 = 0.5 
•when the damping in the absorber is taken as 5 = 0.1. Hie critical 

■value of 6 for second region becomes approximately 0,36 with 5 = 0 . 1 . 

ELgire 4,12 shows the instabili'ty regions of a damped-clamped 
pipe with the absorber attached to it. It is seen that like in other 
cases, with an undamped absorber, two regioas of instability are obtained. 
The damping has more pronounced effect on the instability regions of the 
pipe under cmsideration in comparison to its effect on the instability 
regions of the pipe -with other end ccnditions, lo unstable xegim is 
obtained if the damping in the absorber is takaa as = 0,1 with other 
parameters being same as shown in Hgire 4,12, 











orber attochcd to it.uo 
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4.4.6 Instability Kegicns of a Iteriodically Supported Hpe 

Sie regicsas of instability of a two span pipe, with identical, 
undamped absorbers attached at the mid-point of each bay have beoa 
determined, Ihe wave approach outlined in section 4.3 has bear used 
to comprte the regims of instability. 

It has already bem stated in Chapter 3 that the instability 
regions of a periodically supported two span pipe ccaisist of those of 
a pinned-pinned and a clamped-pinned pipe. Likewise, the instability 
regions for sudi two span pipes with absorbers can also be obtained from 
the results presented above for the pLoned-pimed and the clamped-pinned 
cases, These results have been found to be in exceUoat agreement with 
those obtained by using the wave approach. 

However for number of spans more than two, the results cannot 
be extrapolated from the results of sin^e span pipes. In such cases, 
all the regions of instability can be obtained by using the method 


discussed in secticn 4,3 
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COTGIU SIGNS 

5,1 Conclusions 

Erom the results presented in Chapters 2,3, and 4, the following 
conoLusions are drawn, 

Ihe dynamics of a periodically supported pipe, conyeyirg fl nja 
at constant velocity and pressure, can be conveniently analysed by 
studying the free harmonic wave propa^tion in a similar infinite 
structure, There exist alternate frequeixjy bands of attenuation and 
propagation of free harmonic waves. 

The effect of rotational stiffnesses at 1±ie supports is to shift 
the start of the propagation bands to hi^er frequencies, i)iiereas tbe 
ends of the propagation bands are independent of the rotational stiffness. 
Both the ends of the propa^tion bands shift to lower frequenci^ with 
increasing pressure and velocity of the fluid. 

Free haimonic waves travel with different velocities in the 
positive and the negative directions. This is due to the presence of 
the Coriolis acceleration term in the equation of motion of the pipe, 
Natural frequencies of a periodically supported pipe can be easily 
determined from the curves of propagation constant. Eropagation 
constant values can also be used to determine the natural frequencies 
of a pipe resting on identical supports (with rotational constraints) 
at regular intervals. In both the above cases, the amcxunt of conputation 
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is independent of the mmber of spans of the pipe, Birtheimore, the 
static buckling of periodically supported pipes can also be studied 
from the values of the propagation constants, 

She response of a periodically sujported infinite pipe to a 
convicted harmonic presaire field can be determined using the ireve 
approach, fhe » coincidence’ occurs at a hi^er frequency with increasit^ 
convection velocity of the pressure field. Moreover, the ’coincidence* 
frequency depends on the direction of the convected pressure field, 
i,e, , whether in the same or in the ojposite direction of the fluid 
flow, Ihe damping in the pipe is more effective in cortrolling the 
response than the damping in the supports. 

Ihen the velocity of the fluid is not constant and is having 
a harmonic fluctuation over and above a corstant mean value, the pipe 
experiences parametric instabilities. 

Out of the three methods presented for determination of these 
instabilities, the first method requires prior knowledge of the mode 
shapes of a periodically supported beam, lliis method also beecmes 
unwieldy for large number of spans in the pipe. The number of modes 
to be considered depends on various factors, such as desired accuracy 
of the results, nxmber of instability regions to be determined, e^c. 

The second method does not require ary mode shape approximation 
and the instability regions associated with all the modes are obtained 
from the zeros of a determinant. The order of this determinant;,, however, 
increases drastically with the rumber of spans of the pipe. 
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iEhe third method using the wave approach is most suitable if the 
number of spans are more than two. A simple graphical approadi can be 
used to deteimine all the regions of instability. In this method, the 
amonat of ecmputation is independent of the mmber of spans of the pipe. 

Hie mass ratio parameter (g) does net have s^W significant effect 
on the regions of the instability. This, ne^ecting this poremeter the 
analysis can be made ccnsiderably simpler without any saci^ifice in the 
accuracy* * 

With increase in the pressure and the velocity of the fluid the 
instability regions become wider and are shifted to lower frequencies* 

Damping in the pipe reduces the extent of the instability regions 
and a finite value of the excitation parameter is required to start the 
instability. Damping has more pronounced effect on the instability 
regions associated with the higher modes. Ifysteretic damping is more 
effective, as cempared to viscous damping, in controlling the instability 
regions. 

The instability regions of a pipe are markedly affected by 
attachment of dynamic absorbers. The wave appioadb can also be used 
to determine the instability regions of a periodically supported pipe 
with absorbers attached to each bay. 

In the same range of frequency, the romber of instability regions 
with an absorber are generally more than those without it. However, the 
widths of these regions are much smaller than the original regions. The 
regions can also be effectively controlled by varying parameters like 
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"the mass ratiOj "the ■fainLog ratio ajad the damping ratio of the absoiher’. 
fhG' trends of the effects of these paraaeters siaggest tie eKistence of 
their optisum values. Hence, depending upon specific requir^ients, the 
optimum values of these parameters can be determined, 

5.2 Becommendations for Pufcure Work 

In view of the results obta i ned in the ptresent thesis, following 
recommendations can be made regarding future work in the area. 

(i) Effect of axial tension and contraction of the pipe, neglected 
in the present work, can be accounted for. 

(ii) Regions of combinational resonances for a periodically supported 
pipe can be investigated. 

(iii) The problem of parametric instabilities of pipes with randcsa 
fluctuations in the fluid velocity can be studied. 

(iv) With a suitable choice of an objective ^Unction, optimisation of 
the absorber.^, parameters for the best control of the parametric 
instabilities of pipes can be attempted. 
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APPEEDIX 1 


CAECUEATIOE Op RECEEEAECES POR A PIPE COHWIEO- HaJED 


Consider a pipe of length I , renting on transversely rigid 
supports and let the rotational stiffness at each end be k ^2 
(ligure A-I.l), This foms a periodic element of the infinite pipe- 
line shown in iigure 2,1. The equation for transverse motion of 
the pipe, in absence of any external loading, is [ 57 ] 


El 




(iiu/+ 






+ ]iu 


o » p .y 

8t 3x 


V 


(A— 1* l) 


■sdiere BI is the flezural stiffness of the .pipe. 


and 


y is the displac^ent of Ihe pipe in the transverse direction, 
X is the co-ordinate along the length of the pipe, 

is the mass of the fluid per unit length of the pipe, 

V is the velocity of the fluid, 
p^ is the pressure of the fluid, 

A is the cross sectional area of the pipe, 

Jr 

m is the mass of the pipe per unit length, 

P 


Equation (A-1 ,i) was derived with the following assumptions: 

(i) The pipe behaves like a Euler beam— type structure , 

(ii) The supports are such that thby produce no axial force in the pipe, 

(iii) The axial contraction of the pipe is negligible. 
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Mae significance of various terms in equation (a-1.i) can be 
exia.ained as follows; 

The leading term in equation (A— 1.l) arises due to the heeding 

2 

of the pipe. The term m represents the inertia force of the' 

2 ^ 

pipe and pA is due to the pressure of the fluid. The terms 

2 9x 2 

and 2iii» V are due to the centrifu^l and the 

f g^2 f 3X8t 

Coriolis acceleration of the fluid (moving in a curved iH'th), respectively, 

for 

The acceleration due to ^unsteady flow of the fluid is accountediby 

the term “f » Blnally, the term m^ is the inertia force 

3*fe 

due to the dispLacemexifc of the fluid in the transverse direction* 


If the velocity and pressixre of the fluid are assuined to be conefcanfc 
eveiyvhere *. then tt' = 0* and eqiiation reduces to 

' O w 

El ^ + Cm/t p/p) ^ + anf T + (mf+ V ^ ° 


3x 


3x 


— jyyit 

R)r harmonic solutions, by substituting y = yle in equation 
(A-1,2), one gets the following equation in terms of the ix)&- 
dimensional quantities. 


+ (u^+y) + i 2euQ - o^y = 0 , 






m 


^^ 1/2 


vAiere g = "I , u = (•^) 7SL , y = 


El 


I P 


tdJl 


(A-1.5) 
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Hie solution of equation (A-1,3) can be written as 


4 ix 5 

?(E) = Z ^ “ , 


rt =1 


(A-1.4) 


wiaere the s are the roots of the polynomial 


4 2 2 2 

X - (u + y)x -23unx-fi =0 


{A-1.5) 


To determine the receptances 3^ and 3^^^, used in section 
2*2*1, apply a raon-dimensional unit harmonic momenfc at eM A, ?4ien 
the boundary conditions are 


y(o) 

= 0, 

(A-1,6a) 

y(i) 

= 0, 

(A-1*6b) 

- ?'(o) 

= 1 Kj,y’(o), 

(A-1,6c) 

- ?'(i) 

= ^ y*(i) , 

(A-1,6d) 


^er^ 


h 0 
_ r ^ 

•'r “ El • 


Upon using equation (a- 1,4), equations (A-1,6) yield a system of 
equations 


QG = R , 


where 


Q = 


iX, 


2 <T 

Xi + iX^ 2 


iX. 




1 

iX, 


2 . <T 

^2 ^^2 ~2 


O 


1 

iX. 
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Xj"^ iXj 2 


(a-1.7) 
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X + iX. -4 
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Mow, the reoeptaoces and 


are gtvsn as 


^ = y*(o) 




I 

n=l 

4 iX„ 

I <=n(ig= » 


(4>“1»8a) 


(jUl,8b) 


Diiere liie G^*s are solutions o^ equation (a~i<,7). 

To calculate the receptances and apply a noii- 

diaasnslonal unit harmonic moment at end B. Upon using ■Kie prroper 
boundary conditior^t, the syst^a of equation obtained in this case 
is 


QO = E», 


(A^1.9) 


Ydaere 


fM 

0 

0 \ 


R' = i 


33ie receptances and g^^ are then given as 


BB 


4 


n=1 

''hb = 1 


(A-l.lOa) 


np=1 


» 


(A-l.lOb) 
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Yhere the C^’s are solutioijs of equation (a- 1 . 9 ). 

Star a uniform beam, the receptaaces are related as g , , = - 

AA ijJD 

a-bd ^AB [ 5 ] . Bie to the presence of the Coriolis 

acceleration, this symmetry of the receptaaces no longer holds good 
in case of a pipe. It is found that for a pipe conveying fluid, 

3^ = - and 3 *^ = - 3 ^ , where * denotes the oomilex conjugate. 





AEEEHDIX 2 


DEDEEMIIMIOT OF MTOEAL FHECPEICIES OF A TWO SmN UPE COWErilTG 
FLUIE USme FHEgJESrOT IMJEBMIHAHE 


Gousider a two span pipe resting on equi-spaced supports without 
rotational stiffness as siiown in Figure A-»2*1. Using equation (a~ 1,4), 
the dispLacemert of the pipe, in noo-dimensioml form, can be written 
as 

rl. 


i ^ 1 

y'(cj = ^ c * e ' for span EF, 
rfcl 


-Z(f\ r ^^n^2 

and y (5) = 2. G„ e 


n=1 


n 


for span FG, 


(a^^» la) 


(A^.1b) 


where g and g„ are the non-dimensional co-03rdinates along the 

<0 ^ ^ ■ d. 

length of the pipe for spans EF acd FG-, respectively. The 
12 

coefficients s are unfcmwiK to be determined from boundaiy 

conditions. 

The boundaiy conditions are 


for zero deflection at all the supports 


y\o) = 0 , 

(A-2.2a) 

!l 

O 

(iU2.2b) 

y^(0) = 0 , 

(A-2.2c) 

y^(l) = 0 . 

(A-2.2d) 

for zero momerts at the supports B and G 
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r'(0) =0 , (jU 2 . 2 e) 

y'^(l) =0 . {A-2.2f) 

(c) for coixtijauity of slope and momenb at the support P 

y”(i) = ?2(o) , 

- P'Ho) . <4-2.a) 


Substitution of equation (a— 2« 
a i^stem of homogeneous equations 

Q'C =0 

viiexe 



in equations (a- 2.2) yieldts 

(A-2.3) 



and t denotes the transpose. 
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Por norwtriTial solution of equations (a- 2.3), one gets 

1Q' 1=0 (a*2.4) 

Erequencies satisfying eqimtion (A-2,4) give the mtural 
frequencies of the two span pipe. 

It is seen that for the ease of a two span pipe, the order of 
the deteminarrt |q*I is 8. In gemral, for an N-span pipe this 
order will be 4H. (Uius, the method becomes utsdeldy for large 


values of H. 



APPMDIX 3 


SOME PEOPEEPIES OP THE MATHIEQ-HILL EgjATIOH 

Coiisider the diffeiential eqj&tion 
2 

—I' + fi [ 1 - 26 f(t)] y = 0, (A-3 .i) 

dt 

vhere f(t) is a periodic function with a period 



Since <|) (t + P) =4)(t), equation (A-3.1) does not chan^ its form on 
addition of the period T to t. Therefore, if y(t) is a solution of 
eqaation (A-3 v 1) then y(t + T) is also its solution. 

72 linearly indepeMeafc solutions of 
equation (a-3.i), theny^(t+T) and y^Ct-id!) are also its solutions, and 
consequently can be pEresenfeed in the form of a linear combination a£ 
solution y^ (t) and y^Ct). 

y^(t+T) = a^^ y^(t) + a ^2 ’ (A-3.2a) 

y^Ct+T) ^^22^2^^^" (A-3.2b) 

viiere a. . are constants. 

11 

Thus, the addition of the period T to t results in a linear 
transformation of the initial system of solutions. If one takes some 
other linearly independent solutions instead of the initially cAiosen 
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solutions of y^(t) and YgCt), then the coeffLcdents of the transformtion 
in (A-«3.2), generally, will change. In parti cailar, one can try to choose 
solutions y’^i^t) and y*(t) such that the coefficients a^^ ®’2-j 

equation (1-3.2) vanish. Olie transfonoation in this cas« will take 
its simple^ form and to. 11 be reduced to the nulti plication of initially 
chosen solutions by certain constants'^. 


y*(t-Hl) = y*(t) 


y|(t+!E) = pg y*(t) 


(A-5.3a) 

(A-3.3b) 


Uow, the transformation of the type (A— 3.2) can be reduced to the 
diagonal foim (A-3,3), vhere the constants p.| 2 determined from the 


dmracteristics equation 


^ 

^12 

^21 

C\J 


= 0 


which can be written in the following form 
2 

p _ 2AP + B = 0 
where ^ ^ + ^22^ ’ 


and B = a.,^ - a.,^ ^21' 


(A-3.4) 


Let y.j(t) and y2(t) be two linearly independent solutions of 
equation (a-3.i) satisrfying the initial conditions 


t These are known as Hoquet solutions. 
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y^(o)=a^ , 72 ( 0 )= a^, 

y/o) = 3-, , TgCo) = gg j 


(A^3.5) 


where dots denote differentiation with respect to t. 

Wow, use of equations (A-3.5) in egnations (A-3.2) remits in 


“1 ^2“2 

^1 «1 * ^2®2 
*21 “1 ^ *22 “2 “ 5 ' 2 (®)- 
*21 *1 *22 »2 = J-aCT) 


(A-5,6a) 
(A-3.6b) 
(a-5.6c) 
(a-5 .64) 


and 


Erom equations (a- 3.6), -constants a. . can be detenairaad as 


^t1 = 


“2 “ ^2 

“2 “ “1 ^2 


_ 3., y^(3?) - y^(0}) 

^12 “2 ~ “1 ®2 * 

Op yp(T) - gp yp(T) 

*ai°— s;r ' . - ;eV 

3^ ypCa?) - ypCi) 

"22 " ■“■^Tv^TV 


rvi-c 

Using the values of a. . » constant B can be written as 

f ^ A 

® ~ ^11 ^22 “ ^12 ^21 


or 


B = 


y^(T) ypC^) - ypCT) y|(T) 


a $2 *^2 ^ 1 


(A-5 .7) 



m 

%e 

It can be shown that^constarrfc B is always e(pai to unity. Since 
are solutions of equation (A_ 3 .i), then 

^o f"* ~ = 0 (A-3.8a) 

72 +^0 [ ^ ] ya = ^ (A-.3.8b) 

Multiplying the first identity (A-3,8a) by y 2 ('t), aM the second identity 
by y^(t) , and subtracting one from the other, one gets 

yi('fc) 72^'^) - ° 

on integrating it, one obtains 

y^('fc) y2('*^) ya^'*^) ° > (a-5.9) 

vriaere C is a constant. 

Using initial conditions (at t = o), C can be deterndned as 


0 = .., 

and at t = I, equation (A-3,9) gives 


y^(3?) y2(^) - ya^®) y-j^^) = c = Bg “ “2 


S — ”]# 

Ihus, the characteristic equation (A-3.4) takes the form 


P 


2 


2Ap + 1 = 0 , 


(A^3.10) 


and its roots are related by 
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It has already been shown that among the partimlar soiitions of 
(A-.3.I) two linear independent solutions y* (t) exist whicb satis^ 
equations (a~3.3) 

yJCt+T) = yj(t), k = 1,2. 

■ 5% 

Ihese solutions, vihich acquire a constant maltiilier by the addition of 
the period T to t, can be represented in the form 

, k = 1,2 (a^3.12) 

vdiere ^(t) are certain periodic functions of period !P. 

It follov\e from equation (A-3.12) that the behaviour of the 
soiitions as t " depends on the value of characteristic jcoots (moire 
precisely, on the value of its moduli), la fact, taking iixto accaint that 

la/* = la Jp( + i arg p , 

one can write equation (A-3.12) in the following form 

(t/T)ln}p.| 

y^(-fc) = f^j^(t) e ^ , k = 1,2, (A-5.13) 

where f,, (t) is the bounded function 
1k 

= X^(t) arg p_ 

If the characteristic nunber Pjj. is greater than unity, then the 
corresponding solution (A-3.13) increases unboundedly with time. On the 
other hand, a characteristic mulripLier less than unity implies a decaying 
solution. Finally, if the characteristic number is equal to unity, then 
the solution is periodic, i.e, , it will be bounded in tine. 
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U| > 1, then, as oan be seen from equation (jW 5*10), the 
charaoteristio roots 'wd.ll be real, and one of thoa will be greater than 
unity. In this ease the general solution of equation (a-^.i) will 
unboundedly inci'eafi vdlh time. 

If |l| <1, the characteristic equat;^n has conju^te OOTiplex 
roots and since their jroduct mst be mity, their moduli will be equal 
to unity. The case of ocMiplex characteristic roots corresponds to the 
region of bounded solutions. On the boundaries separating the r^ions 
of bounded solutions, condition |a( = 1 mist be satisfied f 7] . 

Por 111 = 1, multiple roots occur, moreover, it follows from 
equation (a- 3*11), such roots can be either P ^ ^ or P-] “ P2 “ 

In the first case, as seen from equation (A— 3.3), the solution will be 

2 IT 

periodic with a period T = ; whereas in the second case 1iie period 

will be 2T, 

jEierefore, the regions of mbounded solutions are separated from 
the regions of bounded solutions by the periodic elutions witii period T 
anfl 2T, In other words, two solutionsof identical period bound the region 
instability, and two solutions of different periods bound the region of 
stability. 

Now, consider the differential equation of the form 

+ 2e ^ [l-26f(t)] y=0 (a-3.14) 

dt 

where the second teim represents viscous damping on the motion descsrihed 
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Bolotin [ 7 ] tias shown thai the problem of finding the regions 
of instability for equation (a- 3»14) is reduced to the determination 
of the conditions under which it has periodic sblutions with period T 
and 2T, Here also, two solutions of an identical period bound the 
region of unboundedly increasing solutions and tiro solfi^ions of 
different periods bound the region of damped sohitionB. 



APPBJDIX 4 


OP THE MATRICES D , D , D , ME Z 
pp’ cc' cp’ a 

(a) Matrix E ; 

PP 

The elements of the matrix E (order 8 ><8) are 

PP 

^1,3 =’■ 


" ’•13’ 


■5 = ® 

13 




s_. 

+ o ^ A ) e 
4,0 Mo 8 a o^ 

^ /2 

d = S e ^ , 

4,0+4 8P a o » 


<L.=-^r a^B e 
3||3 8 a o 

liA» 

.5 . 1 


N 3 /= 


(L = (X^. + ^ r £5 A ) e 
P»3+4 2o 8 a ^o^ 


20 


^6,3H4 ' ■> ’ 


’‘23’ 

^jd'H ~ ^23 ' 




with 0 = 1»2,3, aad 4. The renaming terms of the matrix E^^ are zero. 

(h) Matrix E : 

' ce 

The elements of the matrix B are same as the elements of the 

cc 

matrix E but for the following changes 
PP 
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“ ^ 13 ’ 

vdth j = 1,2,3, md 4. 

(c) Matrix D ^ ; 

cp 

Ihe elements of the matrix D (order 16 x 16) are 

cp 


d. 1, 


‘^2,3+4 “ 

N/' 


= . 23 


^.72 


"4,3-M2 = ^ 


*^,5 ^io’ 




*^6,^+4 


^,jdS ^id ’ 

2 V" 

^,3+12 = ^23 " 
^,3“° > 

^9,3^6 "■'’ 

X /2 

d =e2j 

10,3-H ® 


*^ 10 , 3+12 


<^4-, • = ^4- ® 

11,3 13 


*13^= 


Si,3-+8 “ ^13’ 
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A „ y2 

d =* A e 

12,0+4 20 » 


d^- . = A , 

12,0+12 2o^ 

2 

d_ . = Af. e 
13,0 10 


S,d-+6 "Nd» 

d = a 2 

‘"l4,d+4 20 


^4,0+12 -^2o’ 


d_ . = (A^ . + y r a) e 

15,0 I0 4 a O'' 

1 2 ^21''*^^ 
d_ . = y r Ci‘^ B e , 

15,0+4 4 a 0 » 




15, d-f6 " “^1d' 

^ ./2 

d , . = - -J- r B e , 

16, g 4 a 0 ^ 




*^16,0+12 “^2o* 


Y/ith 0 = 1,2,3, and 4. OSie remaining terms of the matrix 15 

c 

are zero. 

(d) Matrix : 

Ihe elements of the matrix Z (order 8x8) are 

a 

z . . = 1 , 

1»3 


^2,0 “ ^id’ 



seeci! 




A /2 

^ /2 

z, .= e l/ 


"5,J+4=*'1» 

^/2 

= . 






IP 


"' 6,344 ~ “ 






1j 


"V,j+4 ~ “ 


13' 


- j. 1 2 k ^4^2 

,3 - ‘^3 V , 

3 


^,3+4 ''ij» 

with j = 1,2,3, aad 4. Ihe ronaining teims of the mtriz 


